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SUMMARY 


A  comprehensive  analytical  design  procedure  for  the  installation  of  simple 
pendulums  on  the  blades  of  a  helicopter  rotor  to  suppress  the  root  reactions  is 
presented.  To  achieve  this  goal,  a  frequency  response  analysis  is  conducted  of  typical 
rotor  blades  excited  by  a  harmonic  variation  of  spanwise  airload  distributions  as  well  as 
a  concentrated  load  at  the  tip. 

The  structural  modeling  of  the  blade  includes  elastic  degrees  of  freedom  in  flap  and 
lead-lag  bending  plus  torsion  with  a  hingeless  hub  constraint.  Simple  pendulum 
absorbers  are  individually  considered  for  both  flap  and  lead-lag  types  of  motion.  The 
inertial  reaction  forces  and  moments  at  the  pendulum  hinge  are  formulated.  On  the 
basis  of  a  rational  ordering  scheme  the  general  nonlinear  equations  of  motion  for  the 
rotor-pendulum  system  are  linearized  in  the  perturbation  elastic  blade  displacements 
and  pendulum  angle.  These  linearized  equations  include  the  effects  of  spanwise  airload 
distributions  associated  with  the  elastic  deformations  and  the  cyclic  pitch  environment 
of  forward  flight.  A  quasi-steady  aerodynamic  representation  is  utilized  in  the 
formation  of  these  airloads.  The  solution  of  the  system  equations  is  based  on  their 
representation  as  a  transfer  matrix. 

The  numerical  results  presented  here  pertain  to  both  uniform  and  nonuniform 
hingeless  rotor  blades.  These  results  include  the  effect  of  pendulum  tuning  on  the 
minimization  of  the  hub  reactions.  It  is  found  that  a  properly  designed  flapping 
pendulum  attenuates  the  root  out-of-plane  force  and  moment  whereas  the  optimum 
designed  lead-lag  pendulum  attenuates  the  root  in-plane  reactions.  A  properly  tuned 
pendulum  can  attenuate  the  vibratory  loads  by  generating  appropriate  forces  at  its 
attachment  point  with  the  blade.  These  forces  redistribute  the  loads  on  the  blade  so 
that  only  a  small  portion  of  the  reactions  are  transmitted  to  the  hub. 

For  optimum  pendulum  tuning  the  parameters  to  be  determined  are  the  pendulum 
uncoupled  natural  frequency,  the  pendulum  spanwise  location  and  its  mass.  It  is  found 
that  the  optimum  pendulum  frequency  is  in  the  vicinity  of  the  excitation  frequency.  In 
the  case  of  an  off-tuned  pendulum  the  result  can  be  either  a  slight-to-moderate 
degradation  in  pendulum  absorber  effectiveness  or  a  drastic  increase  in  hub  reactions. 
Although  the  uncoupled  natural  frequency  of  the  pendulum  is  independent  of  the  mass  of 
the  pendulum,  a  proper  choice  of  the  mass  is  required  to  generate  an  optimum  force  to 
attenuate  the  root  reactions. 

A  pendulum  can  be  tuned  and  its  optimum  mass  determined  by  excitation  with  a 
concentrated  simple  harmonic  load  at  the  tip.  However,  it  is  necessary  to  utilize 
distributed  airloads  to  accurately  determine  the  attenuation  of  the  root  reactions.  The 
damping  at  the  pendulum  hinge  has  a  small  effect  on  the  hub  reactions  once  the 
optimum  pendulum  tuning  is  established.  For  the  optimum  pendulum  a  parametric  study 
is  conducted.  The  parameters  varied  include  prepitch,  pretwist,  precone  and  pendulum 
hinge  offset. 


INTRODUCTION 

Vibration  has  always  been  a  significant  detriment  to  helicopters  which  is  not 
surprising  in  view  of  the  fact  that  a  major  percentage  of  a  helicopter's  mass  consists  of 
rotating  components  such  as  engines,  rotors,  gearboxes  and  drives.  These  revolve  at 
different  frequencies,  though  most  are  multiples  or  submultiples  of  the  main  rotor 
frequency,  and  any  out-of-balance  forces  are  easily  transmitted  to  the  light  airframe. 
There  is  also  the  possibility  of  damaging  beat  frequencies  arising  in  those  few  cases 
where  rotational  rates  are  not  synchronous. 


In  addition  to  the  inertially  induced  vibrations,  the  unsteady  air  loads  on  the  main 
and  auxiliary  rotors  significantly  contribute  to  the  vibratory  environment.  These 
unsteady  aerodynamic  reactions  are  caused  by  the  local  periodic  dynamic  pressure  on 
the  main  rotor  blades  in  forward  flight  and  the  cyclic  pitch  changes  imposed  by  the 
rotor  head.  In  response  to  these  unsteady  loads  the  elastic  blades  are  continually 
deforming  which  induces  additional  unsteady  aerodynamic  loading. 

An  increasing  demand  for  the  reduction  of  helicopter  vibration  has  been  dictated  by 
the  expansion  of  flight  envelopes  coupled  with  more  stringent  requirements  for  crew 
and  passenger  comfort  as  well  as  improved  reliability  and  maintainability. 

The  rotor  system,  which  transmits  the  vibratory  airloads  to  the  fuselage  through 
the  rotor  shaft,  is  one  of  the  most  significant  contributors  to  the  vibrations  of  the 
helicopter.  The  rotating  blades  create  vibratory  airloads  containing  all  harmonics  of 
the  rotor  rotational  frequency.  As  the  loads  from  the  individual  blades  combine  at  the 
rotor  hub,  some  harmonics  cancel  each  other  while  others  are  additive.  The  loads  that 
are  additive  are  passed  from  the  blades  to  the  pylon  and  then  to  the  cabin  through 
complicated  dynamic  load  paths.  These  loads  are  felt  as  vibratory  forces  and  moments 
whose  frequencies  are  integer  multiples  of  the  blade  passage  frequency  (number  of 
blades  times  rotational  frequency). 

Dynamic  systems  that  reduce  the  vibration  in  helicopters  may  be  classified  into 
five  groups: 

(1)  Excitation  reducers  (blade  aerodynamics) 

(2)  Attenuators  (blade  dynamics  with  or 

without  hub  motion) 

(3)  Absorbers  (hub  dynamics) 

W  Isolators  (pylon  dynamics) 

(5)  Cabin  Suppressors  (fuselage  dynamics) 

These  systems  are  described  below  and  specific  examples  are  given. 

Excitation  reducers  are  devices  that  either  reduce  aerodynamic  forces  on  the 
blades,  or  alter  the  hub  moments  and  shears  by  generating  counter-influencing 
aerodynamic  forces  on  the  blades.  Higher  harmonic  blade  pitch  control  devices  (ref.  1) 
fall  into  this  category.  In  reference  1  a  wind  tunnel  test  has  been  used  to  evaluate 
higher  harmonic  blade  pitch  for  the  reduction  of  helicopter  vibration.  The  higher 
harmonic  pitch  was  obtained  by  higher  harmonic  oscillation  of  the  swashplate.  The  test 
results  indicated  reduction  of  fuselage  vibratory  loads. 

The  attenuators  are  those  dynamic  devices  that  result  in  low  hub  shears  and 
moments  while  the  blades  are  subjected  to  specified  aerodynamic  excitation  forces.  An 
example  of  this  category  is  a  midspan  or  tip  weight  that  may  move  the  natural 
frequency  of  a  mode  away  from  exciting  harmonics  or  reduce  the  modal  forces  at  the 
blade  root. 

Absorbers  are  devices  that  are  incorporated  at  the  hub  (or  blade  root)  and  absorb  a 
large  portion  of  the  hub  shears  or  moments  causing  transmittal  of  only  the  remaining 
forces  to  the  mast.  Examples  are  the  simple  and  bifiiar  pendulums. 

Isolators  are  generally  integral  to  die  pylon  assembly  and  reduce  greatly  the  forces 
and  moments  transmitted  to  the  fuselage  with  respect  to  the  forces  applied  on  top  of 
the  mast. 

Cabin  suppressors  reduce  the  vibration  in  the  cabin  for  prescribed  forces  and 
moments  applied  by  the  pylon  on  the  fuselage. 

One  of  the  five  categories,  the  pendulum  absorber,  has  proven  to  be  quite 
successful  as  reported  in  references  2  to  6.  Such  a  technique  consists  of  mounting  a 
simple  or  bifiiar  pendulum  on  each  rotor  blade  near  the  hub.  The  blade-mounted 
pendulum  absorbers  modify  the  response  dynamics  of  the  blades,  and  reduce  the 
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helicopter  vibrational  levels.  All  reported  work  has  consisted  of  an  expensive  and  time- 
consuming  flight  test  program  to  establish  a  set  of  acceptable  pendulum  parameters. 
As  the  flight  tests  become  more  and  more  expensive,  an  analytical  study  is  needed  for 
pendulum  design. 

Murthy  and  Hammond  (ref.  7)  have  studied  the  effect  of  a  pendulum  absorber  on 
the  natural  frequencies  and  mode  shapes  of  a  rotor  blade.  A  vibration  analysis  of  rotor 
blades  undergoing  coupled  flapwise  bending,  chordwise  bending,  and  torsional  vibrations 
with  a  spherical  pendulum  absorber  was  presented.  A  portion  of  the  investigation 
included  a  reduction  of  the  spherical  pendulum  to  a  simple  flapping  pendulum.  It  was 
shown  that  the  pendulum  absorber  displaces  original  blade  natural  frequencies  that  are 
in  the  vicinity  of  the  pendulum  tuning  frequency(ies).  The  spherical  pendulum 
introduces  two  modes  between  the  displaced  original  modes  of  the  blade,  whereas  the 
simple  flapping  pendulum  introduces  one  mode.  Mode  shape  results  were  not  presented 
in  the  report.  It  is  worth  mentioning  that  the  equations  of  motion  were  based  on  the 
linear  analysis  presented  by  Houbolt  and  Brooks  (ref.  8),  and  no  damping  terms  were 
included.  Also  the  damping  at  the  pendulum  hinge  was  not  considered. 

The  principal  aim  of  the  present  study  is  to  establish  comprehensive  analytical 
design  procedures  and  criteria  for  the  installation  of  simple  pendulums  on  the  blades  of 
a  helicopter  rotor  to  suppress  the  force  and  moment  reactions  at  the  hub  during  forward 
flight  conditions.  ,7,-r 

-The  objectives  of  the  present  investigation  can  be  stated  as  follow^: 

(1)  Develop  a  mathematical  model  to  represent  the  blade-pendulum  system.  A 
single  nonuniform  rotor  blade  with  a  hingeless  hub  restraint  undergoing  coupled  flapwise 
bending,  chordwise  bending,  and  torsional  vibrations  is  considered.  Simple  pendulum 
absorbers  are  individually  treated  for  both  flap  and  lead-lag  types  of  motion.  The  blade 
is  excited  by  an  azimuthal  harmonic  variation  of  spanwise  airload  distributions 
associated  with  the  elastic  deformations  and  the  cyclic  pitch  environment  of  forward 
j<ight. 

(2)  Find  the  dynamic  response  characteristics  of  the  blade-pendulum  system,  using 
the  transfer  matrix  method  (ref.  9). 

(3)  Determine  the  optimum  pendulum  tuning  to  suppress  the  hub  reactions.  This 
entails  the  minimization  of  these  reactions  by  appropriate  variations  of  the  pendulum 
parameters  for  a  given  excitation  frequency.  The  pendulum  parameters  include  the 
uncoupled  pendulum  frequency,  hinge  spanwise  loaction  and  pendulum  mass. 

(4)  Conduct  a  parametric  study  of  the  optimum  tuned  configuration.  The 
parameters  to  be  varied  include  pendulum  hinge  offset,  precone,  prepitch  and  pretwist. 
The  intention  of  this  investigation  is  to  document  the  effects  of  these  parameters  on 
the  optimum  configuration  previously  established  and  thereby  provide  useful  design 
criteria  for  future  installations  of  pendulum  absorbers.  v 
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SYMBOLS 

a  airfoil  lift-curve  slope 

60  reference  frame  which  rotates  with  speed  Q  with 

respect  to  the  stationary  inertial  frame 

C  pendulum  hinge  friction  coefficient,  equation  (35) 

Cj  rotor  thrust  coefficient 

c  blade  chord 
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airfoil  profile  drag  coefficient 
Young's  modulus  of  elasticity 

mass  centroid  offset  from  elastic  axis,  positive  when 

forward  of  the  elastic  axis 

aerodynamic  loads  vector,  equation  (133) 

components  of  aerodynamic  forces  per  unit  length  in  y', 

y',  z'  directions 

shear  modulus 

acceleration  due  to  gravity 

"fixed”  unit  vectors  with  &  in  direction  of  rotation  of 
figure  A-l 

blade  cross-section  moment  of  inertia  about  y'  axis 
blade  cross-section  moment  of  inertia  about  z'  axis 
unit  vectors  assoicated  with  undeformed  blade  coordinate 
system,  figure  A-l 

unit  vectors  associated  with  deformed  blade  coordinate 

system,  figures  1,  A-2 

torsional  rigidity  instant 

radius  of  gyration  of  blade  cross  section 

mass  radius  of  g;  ution  of  blade  cross  section 

mass  radii  of  gyration  about  n  and  t  axis,  respectively 

components  of  generalized  nonconservative  (aerodynamic) 

forces  in  y  and  z  directions 

length  of  pendulum 

mass  of  the  pendulum 

twisting  moment  about  the  x  and  x*  axes 

bending  moment  about  y  and  y'  axes 

bending  moment  about  z  and  z'  axes 

generalized  nonconservative  (aerodynamic)  moment  about 

x  axis 

mass  per  unit  length  of  the  blade 

time  independent  average  normal  force  on  the  pendulum 

bearing,  equation  (33) 

fixed  point  irt^and^fat  the  root  of  the  blade,  figure  A-l 
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point  transfer  matrix  of  the  pendulum 

length  of  blade 

inertial  reference  frame 

transfer  matrix  of  the  blade  from  0  to  x 

initial  radial  blade  tension,  equation  (10) 

time 

radial,  tangential,  and  perpendicular  components  of 
velocity  for  blade  airfoil  section 
elastic  displacements  in  the  x,y,z  directions, 
respectively,  figure  A-2 
wind  velocity 

axial  force  (tension)  in  the  x,  x*  directions, 
respectively 

shear  force  in  the  y,  y'  directions,  respectively 
shear  force  in  the  z,  z'  directions,  respectively 
induced  downwash  velocity  at  rotor,  positive 
downward 
pendulum  weight 

inertial  axis  system  with  origin  at  hub  centerline  and 
Z  normal  to  hub  plane,  figure  A-l 
mutually  perpendicular  axis  system  with  x  along  the 
undeformed  blade  and  y  toward  the  leading  edge, 
figures  A-l 

a  set  of  Cartesian  coordinates  with  origin  at  a  point 
along  the  deformed  blade  with  x'  remaining  tangent 
to  the  elastic  axis  and  y'  and  z'  along  principal  axes 
for  the  cross  section,  figures  1,  A-2 
state  vector,  equation  (1) 
shaft  angle 

precone  angle,  figure  A-l 

small  parameter  of  the  order  of  the  bending  slopes, 
also  airfoil  section  pitch  angle  with  respect  to  free- 
stream  velocity,  figure  7 

sectional  coordinate  normal  to  n  axis  at  elastic  axis 
(same  as  z'),  figure  1,  also  damping  ratio 
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sectional  coordinate  corresponding  to  major  principal 
axis  for  a  given  point  on  the  elastic  axis  (same  as  y’), 
figure  1 

pitch  angle,  equation  (125) 
collective  pitch,  equation(i25) 
cyclic  pitch  components 
pretwist,  equation  (126) 

pretwist  per  unit  span,  positive  when  tip  angle  is 
larger,  equation  (126) 
inflow  ratio,  positive  upward 

slope  of  deflection  curve  in  the  plane  of  rotation, 
dv/dx 

advance  ratio,  OR  =  Vj  cos  a 

slope  of  deflection  curve  normal  to  the  plane  of 

rotation,  dw/dx 

mass  density  of  blade,  also  mass  density  of  air 

elastic  twist  about  the  elastic  axis 

blade  azimuth  angle  measured  from  downwind  position 

in  direction  of  rotation 

rotational  speed  of  the  blade 

uncoupled  natural  frequency  of  the  simple  pendulum 


at  the  pendulum  hinge  location 
at  the  pendulum  mass  location 

immediate  left 
immediate  right 
transpose  of  matrix 
inverse  of  matrix 

differentiation  with  respect  to  x,  also  deformed 
coordinate  system 
differentiation  with  respect  to  time 
amplitude  of  simple  harmonic  quantity 
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BLADE  EQUATIONS  OF  MOTION 


The  equations  of  motion  which  are  used  to  represent  the  structural  dynamic 
characteristics  of  the  helicopter  blade  are  based  on  the  development  reported  by 
Hodges  and  Dowell  in  reference  10.  The  mathematical  model  chosen  to  represent  the 
rotor  blade  consists  of  a  straight,  slender,  variably  twisted,  nonuniform  beam  which  can 
undergo  combined  flapwise  bending,  lead-lag  bending  and  torsion.  The  elastic  axis,  the 
mass  axis,  and  the  tension  axis  (area-centroid  axis)  are  taken  to  be  noncoincident.  The 
elastic  axis  is  inclined  to  the  plane  of  rotation  at  a  small  angle  in  order  to  accommodate 
any  built-in  coning  (precone).  The  use  of  an  actual  helicopter  rotor  blade  as  a  model 
would  be  a  more  formidable  task  because  of  the  nonhomogeneous,  anisotropic 
structures  found  in  typical  blades. 

Several  orthogonal  coordinate  systems  and  coordinate  transformations  (Appendix  A) 
are  employed  in  the  derivation  of  the  equations  of  motion.  The  beam  cross  section  is 
shown  in  figure  1  before  and  after  deformation.  The  n  and  £  axes  are  the  principal  axes 
of  the  cross  section  for  the  shear  center.  The  cross  section  is  assumed  to  be  symmetric 
with  respect  to  the  n  axis.  The  deformed  beam  is  shown  in  figure  2(a)  with  force  and 
moment  resultants  acting  on  the  face  of  a  cross  section.  The  x,  y,  z  coordinate  system 
is  the  preconed  undeformed  system.  At  any  point  along  the  deformed  beam,  x'  is 
tangent  to  the  deformed  elastic  axis.  The  y'  and  z'  are  identical  to  the  n  and  C  axes, 
respectively,  when  the  beam  is  deformed. 

A  mathematical  ordering  scheme  which  is  consistent  with  the  assumption  of  a 
slender  beam  is  adopted  for  the  purpose  of  systematically  discarding  elastic  and 
dynamic  terms  which  are  of  higher  order  in  the  equations  of  motion.  The  ordering 
scheme  is  consistent  with  the  small  deformation  approximation  in  which  the  strains  are 
negligible  compared  to  unity  and  the  dimensionless  axial  deflection  u/R  is  generally 
taken  to  be  of  the  same  order  of  magnitude  as  the  square  of  v/R  or  w/R  and  thus  is 
small  with  respect  to  unity.  These  assumptions  can  be  systemized  by  introducing  e ,  a 
parameter  of  order  v/R  or  w/R.  The  ordering  scheme  associated  with  the  above 
consideration  is  summarized  in  Table  1.  Within  the  force  and  moment  expressions, 
terms  of  order  e  l  are  neglected  with  respect  to  ur^ty.  Thus,  if  the  largest  terms  of  the 
force  expression  are  (X  e  j,  then  all  terms  of  CX  e  )  are  retained  (first-order  terms),  all 
terms  of  0(  e;  are  retained  (second-order  terms),  and  generally  terms  of  CX  eb)  are 
discarded. 

The  blade  equations  of  motion  will  be  written  in  terms  of  the  state  vector  {  Z } 
where 

{Z(x,t)}T  =  Lvw*pvMx  My  Mz  Vx  Vy  V2J  ( 1 ) 

It  may  be  noted  that  the  above  state  vector  does  not  include  the  axial  deflection,  u.  In 
order  to  obtain  a  set  of  linear  equations  for  the  state  vector  components,  u  must  be 
eliminated.  The  components  of  the  state  vector,  { z},  can  be  chosen  several  ways,  but 
they  are  chosen  here  such  that  they  represent  the  physical  quantities  of  deflections, 
slopes,  moments,  and  shears.  This  is  not  absolutely  required  but  highly  preferable  for 
application  of  the  transfer  matrix  to  obtain  the  dynamic  response  characteristics. 

The  equilibrium  of  the  forces  and  moments  that  act  on  a  differential  beam  element 
is  now  considered.  In  this  consideration  the  element  is  formed  by  slices  parallel  to  the 
yz-plane,  because  this  choice  leads  to  rather  simple  results.  The  forces  that  act  on  such 
an  element  are  shown  in  figure  2(a);  the  moments,  in  figure  2(b).  The  quantities  p  ,  p  , 
p  ,  q |  ,  a  and  q  are  resultant  force  and  moment  loadings.  Summation  of  the  moments 
aoour  the  x-,  y-,  and  z-directions  lead  to  the  following  equilibrium  conditions  for 
moment  and  shear: 
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M'x  +  Vz  v'  -V  w'  +  qx 

=  0 

(  2a) 

M'  -V  +  V  w'  +  q 

y  z  x  y 

=  0 

(  2b) 

M'z  +  Vy  -vx  v’ +  q2 

=  0 

(  2c) 

V'  +  p 

X  Kx 

=  0 

(  2d) 

V'  +  p 

y  y 

=  0 

(  2e) 

V'  +  p 
z  Kz 

=  0 

(  2  f  ) 

The  applied  loads  p  ,  p  ,  p  ,  ^ue  to  inert*a*  an<^  aerodynamic  loadings  are 

functions  of  v,  w,  <j>  ancf a  l&ngitudinarradial  deformation  u  in  the  x  direction.  The  only 
effective  aerodynamic  loading  to  be  considered  in  this  analysis  will  be  forces  in  the  y 
and  z  directions  given  by  L  and  L  ,  respectively,  and  moment  about  the  x  -  axis  given 
by  M  .  .  The  inertia  loadings,  whi^h  consist  of  the  forces  and  moments  which  oppose 
accelerations  of  the  blade  elements,  are  derived  in  reference  10  whereas  the 
aerodynamic  loadings  L  ,  L^,  and  M.  are  left  in  this  symbolic  form.  These 
aerodynamic  loads  are  derived  later  in  this  text.  The  procedure  that  follows  is 
employed  in  the  derivation  of  the  inertia  loading  which  is  given  in  reference  10.  The 
acceleration  of  any  mass  particle  on  the  vibrating,  rotating  blade  is  derived  and  the 
components  in  the  x-,  y-  and  z-  directions  are  obtained.  These  component  accelerations 
include  terms  for  the  Coriolis,  transverse  and  centripetal  accelerations  along  with 
higher  order  terms.  The  inertia  force  and  moment  loadings  are  then  obtained  by 
integrating  over  the  cross  section.  By  dropping  higher  order  terms  in  accordance  with 
the  adopted  ordering  scheme,  the  following  resultant  loadings  to  second  order  are 
obtained: 


Px 

2 

=  m(  ft  x  +  2SlV)-mBpc'w 

(  3a) 

py 

=  -mV  +  me  $  sin  9  +  m  fl  ^  Ev  +  e  cos  (6  +  4>  0 

(  3b) 

+  2m  fl(  8  w-fl)  +  2m  fie  (v1  cos  0  +  w'sin  0 )  +  Ly 

pz 

....  2 

=  -mw  -me  $cos  0  -m  fl  8pcx  -2mfl8pcv  +  Lw 

(  3c) 

% 

2  2 

=  me(v  sin  0  -w  cos  0 )  -m  fl  ev  sin  0  -m  fl  e  0p(;x  cos  0 

(  3d) 

-mk?£  -A  m  fl^kjl  -U?  )  sin  2(0  +  0 )  +  M  . 
m  i  m-.  m  j  _  o 

qy 

2 

=  me(  fl  x  sin(  0  +  $)  +  2  flv  sin  0  ) 

(  3  e  ) 

^z 

2 

=  -me(  fl  x  cos(0  +  $ )  +  2  flv  cos  0  ) 

(  3f  ) 

The  underlined  higher-order  terms  could  be  dropped  consistent  with  the  ordering 
scheme,  but  this  would  eliminate  torsional  inertia.  In  the  special  case  of  rotor  blade 
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configurations  of  very  low  torsional  rigidity,  these  terms  may  contribute  substantially 
to  the  magnitude  of  the  torsional  natural  frequency;  hence,  they  will  be  retained  (ref. 
10). 

There  are  now  6  equilibrium  equations,  equations  (2a)  -  (2f),  and  10  unknowns,  u,  v, 
w,  4>,M,M,M,V,V,  and  V  .  In  addition,  4  force  and  moment  deformation 
relations  are  required  to  express  M  ,z  M  ,  M  ,  V  in  terms  of  u,  v,  w,  .  The  force  and 
moment  resultants  are  expressed*  in  yterms  oi  the  deformations  by  resolving  the 
distributed  stresses  into  a  resultant  force  and  moment  system  acting  at  the  elastic  axis. 
These  relations  given  in  reference  10  are  written  in  the  deformed  body  axis  system  as 


Vx’  = 


EA  |  u'  +  ^  +  k \  ®  '  ^,-e  a  Lv"  cos  ( ®  +  ♦ )  +  w"  sin  (0  +  <t>  )  J  |  (4a) 


(2  (2 

Mx,  =  GO  <j>'  +  EAk2(0  ♦  $)'  (u*  +■  y  ♦  ~ )  +  EB*j0  'V 


(  4b) 


-EB|  0'(v"  cos  0  +  w"  sin  0 ) 


My,  =  Ely,  j\"  sin  (  0  +  <J) )  -w”  '  os  (0  +  4>)]  (  4c  ) 


Mz,  =  EIZ,  [v"  cos  (0  +  q>)  +  w”  sin  (0  +  <fr)]-EAeA  (u*  +  ^ )  -EB£ 6'  (4d) 


where  and  B*  are  section  constants  defined  by 


Bf  =  JJ (  n2  +  cVdndc 
A 

B2  =yjo(  n2  +  c2)  dn  dc 

A 

and  A  is  the  blade  cross-section  area  effective  in  carrying  tension. 

The  underlined  terms  may  be  neglected  according  to  the  ordering  scheme. 
However,  as  discussed  before,  these  terms  are  important  for  configurations  with  low 
torsion  stiffness.  In  equations  (4a)  -  (4d),  the  terms  which  come  from  the  warp  effect 
are  omitted  due  to  its  relative  smallness  for  closed  cross  sections. 

The  force  and  moment  deformation  relations,  equations  (4a)  -  (4d),  can  be  readily 
transformed  into  the  undeformed  body  axis  system  using  the  transformation  matrix  of 
Appendix  B.  The  resultant  moments  ares 
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Mx  =  G3  +  EAk^  (Q  +  *)'  (u’  +  y  +  ) 

+  EBf  0'2  $'-EB*  0'(v"cos  0  +  w”sin  0) 

2  ,2 

-  EAeA  (u '  +  (v'  sin  0  -w'  cos  0  ) 

+  (EIZ,  -Ely,)  [|  (v'v"  -w'w")  sin  2  0  +  (v V* 

♦  w*v")  sin2  0  1  +  El  ,  v  V’  -El  ,  v"w’ 

y  z 

2  2 

My  =  GJ  4»'v'  +  EBJ  0 ' sin  0  +  EAe^  (u1  +  y-  +  ^-) 

•  sin  (0  +  $)  -(El2,  -EIy()  [w"  sin2  (0  +  $) 

+  \  v"  sin  2(0  +  $)]  -Ely,  w" 

Mz  =  G3  $*w'  -EB£  0’  $'  cos  0  -EAe^u' 

2  2 

♦  ^  +  y- )  cos  (  0  +  $ )  +  (EIz,  -Ely,) 

•  j  [w"  sin  2  (  0  +  $ )  -v"  sin2(0  +  $)]  +  El z,v " 

The  resultant  axial  force  in  the  undeformed  body  axis  system  is 

Vx  =  Vx,  -Vy,  [v'  cos  (  0  +  ♦)  +  w'  sin  ( 6  +♦ )]  -Vz,  [-v’  sin  (0  +$  )  +  w'  cos  (0  + 


(  5a) 


(  5b) 


(  5c) 


♦  >] 
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For  rotor  blade  applications  the  second  and  third  terms  on  the  right-hand  side  are 
negligible  to  second  order  (ref.  10);  thus, 


2  2 

=  Vx,  =  EA  u'  +  Y  +  Y  +  kA  0  '  ♦'  -eA  [v"cos  (6  +  $)  +  w"  sin  (0  +  ♦)]  (  6 ) 

■ 


In  the  above  expressions  all  terms  above  second  order  have  been  neglected.  As  a  final 
step  the  following  definitions  are  introduced. 


w'  =  v  (7  b) 

Equations  (2),  (5),  (6),  and  (7)  represent  12  independent  equations  for  the  11  state 
variables  and  the  elastic  elongation,  u.  In  order  to  obtain  a  set  of  linear  equations  for 
the  state  vector  components,  u  must  be  eliminated.  This  will  be  accomplished  by 
defining  the  time  independent  tension,  T  ,  as 

R 

T  =/  mil2xdx 

°  l 

In  the  process  of  linearization,  the  following  approximations  will  be  made 


sin  (0  +  $)  a  sin  0  +  cos  0 

cos  ( 0  +  $ )  a  cos  0  -  $  sin  0 

j  sin  2  (  0  $)  a  j  sin  2  0  +  $  cos  2  0 

2  2 
sin  (0  *  $)  a  sin  0  +  $  sin  2  0 


(  8) 


Equation  (6)  can  be  rewritten  as 


v«2  v2  Vx  2 

u’  +  -5-  *  —  =  -  kA  0  '  +  eA  y'  cos  0  +  eAv  '  sin  0 


(  9) 


The  left-hand  side  of  equation  (9)  appears  in  equations  (5)  with  both  constant  and  linear 
coefficients.  When  the  coefficient  is  constant  it  shall  be  represented  as  written  above 
but  with  a  linear  coefficient  it  will  be  approximated  as 


,  ,  V2  wl2.  ^  To 
(u  +  2  *  2  *o  “  EA 


(  10) 
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In  the  process  of  linearization,  all  steady-state  values  are  neglected  except  the 
initial  radial  tension,  and  only  the  linear  terms  in  the  state  variables  are  retained. 
These  steady-state  values  would  correspond  to  an  equilibrium  flight  condition  in  hover. 
It  is  anticipated  that  these  equilibrium  hover  values  would  have  little  effect  on  the 
dynamic  response  characteristics  in  the  forward  flight  condition  of  this  study.  The  set 
of  linearized  first-order  differential  equations  are  given  below 


v'  =  y 

( 

11) 

W1  =  V 

( 

12) 

*'  (GJ  +  k 2  T 

A  C 

(  +  EB*  8  ,2)  =  -k2A  6 '  Vx  ♦  EA(k2  6  ')2  ♦' 

-EAkA  0  'eA  (  y'  cos  0  +  v'sin  0 )  +  EB£  0  •(  | 

y' cos  0 

+  v'sin  0)  ♦  eA  TQ  (ysin0  -vcos  0)  +  Mx 

( 

13) 

Vi '  (EIZ,  cos20 

+  Ely,  sin20  )  =  Mz  ♦  EBJ  0  '  cos  0  +  eA 

COS0  Vx 

-  EAkA  eA  0  ' cos  0  +  EAeA  (  y'cos2  0  + 

v  1  sin  0 

cos  0) 

-eA  Tq$  sin  0  -j(EI2,  -EIy^v  '  sin  20 

( 

14) 

V '  (Ely,  cos20 

♦  EIz,  sin20  )  =  -My  +  EB|  0  '  sin0  +  eA  Vx  sin  0 

-EAkA  eA0  *  $'  sin  0  +  EAeA  (  y'sin  0  cos  0 

+  v '  sin 

2e> 

♦  e  A  Tq  4>  cos  0  -  j  (EIZ,  -Ely,)  y  '  sin  2  0 

( 

15) 

2 

M'x  =  m a  e  x(vcos  0  -  y sin  0  )  +  (eA  Tq  cos 6 )'v  -(eA  T 

9  sin0 )' 

M 

+  mk^  $  +me  [(fl2v  -v)  sin  0  +  w  cos  0] 

+  mfi2(k2  -k2  )$  cos  2 0  -M  . 

nrij  m  j  $ 

( 

16) 

M’v  =  V2  -To 

2 

v  -2  mQev  sin  0-mfl  ex  ♦  cos  8 

( 

17) 

13 


(  18) 


2 

=  Tq  y-Vy  +  2m flev  cos  8  -mft  ex  $  sine 
=  -2  m  flv  (19) 


"V 


V 

z 


2  2 

mV  -me  $  sin  0  -m  ft  v  +  mfte$sin0  -2  m  Qg 


pc 


-2  m  fte  ycos  8  -2  mft  e  v  sin  8  -L 


m#  +me  if  cos  0  +  2  m  00  v  -L 

pc  w 


w 


(  20) 

(21) 


It  may  be  noted  that  the  three  equations  lor  the  elastic  slopes  (eqs.  (13),  (14),  and 
(15))  are  coupled  in  their  derivatives  and  must  be  solved  simultaneously.  These 
equations  are  decoupled  by  a  straight  forward  application  of  Cramer's  rule  to  obtain  a 
set  of  linear  first  order  differential  equations  in  terms  of  the  state  vector  so  that  the 
transfer  matrix  method  can  be  applied  as  discussed  later  in  this  study. 

Excitation  of  the  blade  is  due  to  various  simple  harmonic  airloads.  The  solution  to 
the  set  of  equations  will  be  treated  as  individual  contributions  associated  with  each 
harmonic  of  the  excitation  loading.  Such  a  solution  for  the  nth  harmonic  of  the  state 
vector  can  be  written  in  complex  form  as 


{Zn(x,t)}  =  Re({Zn(x)}  eiut)  (22) 

where  {  Zn(x)}  is  the  amplitude  vector  of  the  nth  simple  harmonic  response,  and 
a)  =  nQ  : the  frequency  of  the  excitation  loading 

i  =vr 

The  above  solution  will  be  discussed  in  more  detail  with  respect  to  the  aerodynamic 
loading.  Substituting  equation  (22)  into  equations  (1 1)  to  (21)  yields  the  following  matrix 
equation  for  {  Z] ,  the  nth  harmonic, 


{f(x)}«  =  [  A  (x)]  {  f(x)}  ♦  (Pa(x)} 


where  { Fj  is  the  amplitude  vector  of  the  aerodynamic  loads,  defined  by 


(  23) 


{F,)T  ,[0  0  0  0  0  -rtf,  0  0  0  -Cv  -Lwj 
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Note  that  the  matrix,  [a],  includes  complex  terms  due  to  retaining  the  linear  rate 
terms  (Coriolis  intertial  forces). 

Evaluation  of  the  order  of  magnitude  of  the  individual  elements  of  the  complex 
matrix  £a]  is  made.  Consistent  with  the  adopted  ordering  scheme,  only  terms  of  first 
and  second  order  are  retained.  It  may  be  noted  that  in  the  process  of  decoupling  eqs. 
(13),  (14),  and  (13),  terms  of  higher  order  are  encountered.  The  resulting  nonzero 
elements  of  £a^  become 


A14 

2 

1 

A23 

= 

1 

A34 

= 

EIy,(EIz,-EAe^)  eA  TQ  sin  6  /D 

A35 

= 

-a34  cot  Q 

A36 

= 

EIyI(EI2,-EAeA)/D 

A37 

= 

Eiyl  e  *(EAkA  eA  -EB|P sin  9^D 

eo 

< 

= 

-A yj  COt  e 

A39 

= 

Ely,  9  KEB*  eA  -EIZ,  k^)/D 

A43 

= 

-GJ  eA  Tq  (EIz,  -EA  eA)  sin  0/D 

A46 

= 

A38 

A47 

= 

GJ  (El  .  -El  ,  -EA  e  f)  sin  0  cos  0/D 
z  y  a 

oo 

< 

= 

-A^j  sin  0  /(eA  TQ)  +  GJ  EIyt  cos2  0  /D 

A49 

= 

GJ  El  .  e .  cos0  /D 
y  a 

A53 

s 

-A^3  cot  0 

A56 

= 

-A37 

A37 

= 

2 

A^3  cos  0  cot  0/(eA  T0)  -GJ  Ely,  sin  0 

A38 

= 

"A47 

A39 

= 

\9  tan  0 
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A61 

= 

2  2 

me  ( ft  +  u  )  sin  0 

A62 

= 

2 

-me  a)  cos  6 

A63 

= 

-mk2  u  2  +  m  fl2  (k^  -  k^V  )  cos  2  0 

m  m2  nij 

A60 

= 

2 

-m  fl  ex  sin  0  -(e,  T  sin  0 )' 

A  O 

A65 

= 

2 

m  fl  ex  cos  0  +  (e^  Tq  cos  0  )' 

A71 

= 

-2  im  flue  sin  0 

A73 

= 

2 

-m  fl  ex  cos  0 

A75 

= 

‘To 

A7,ll 

= 

I 

A81 

= 

-A71  cot  0 

00 

< 

a 

2 

-m  fl  ex  sin  0 

A84 

= 

To 

A8,10 

= 

-1 

A91 

= 

-2  im  flu 

A92 

= 

-mu20pc 

A10,l 

= 

2  2v 

-m  ( fl  +  u  ) 

A10,2 

= 

-2  im  fiwSpj. 

AI0,3 

= 

A61 

A10,4 

= 

"A81 

•TN 

•k 

O 

< 

a 

A71 

All,l 

= 

"A10,2 

-mo) 


2 


A 


11,3  ' 


where 


D  =  G 3  EIy,  (EIZ,  -EA  e  2) 


PENDULUM  EQUATIONS 

The  pendulum  configurations  considered  are  of  the  simple  type  which  are  free  to 
oscillate  in  the  flap  and  lead-lag  modes.  The  free  vibration  equations  of  motion  for 
simple  pendulums  in  both  die  flap  and  iead-lag  degrees  of  freedom  are  developed  in  an 
identical  fashion.  These  equations  are  formulated  by  first  determining  the  acceleration 
of  the  pendulum  mass  in  the  inertial  system  and  then  transforming  it  to  the  deformed 
blade  system.  The  absolute  acceleration  is  used  to  express  the  inertial  reaction  force 
and  moment  at  the  pendulum  hinge.  Setting  the  reaction  moment  about  the  pendulum 
hinge  to  zero  yields  the  equation  of  motion.  These  nonlinear  equations  of  motion  are 
subsequently  linearized  by  assuming  small  oscillations  for  the  pendulum  about  the 
steady-state  condition.  By  assuming  a  simple  harmonic  solution  for  the  linearized 
equation  of  motion  of  the  free  pendulum,  the  deflection  of  the  pendulum  is  obtained  in 
terms  of  the  blade  deflection.  The  point  transfer  matrix  at  the  pendulum  hinge  station 
on  the  rotor  blade  is  derived  by  simple  equilibrium  considerations  of  forces  and 
moments  across  the  pendulum  after  linearizing  the  inertial  forces  and  moments. 

Acceleration  Components 

The  acceleration  of  the  pendulum  mass  is  determined  with  respect  to  an  inertial 
reference  system.  This  entails  appropriate  transformations  from  the  inertial  to  the 
rotating  system,  to  the  preconed  system,  and  finally  to  the  deformed  blade  section.  The 
derivation  of  the  acceleration  components  is  given  in  Appendix  B.  Equations  (B-18)  to 
(B-20)  for  these  acceleration  components  are  reproduced  here  as 


ax,  =  5  -2(1  v  -  fl  (x  ♦  £  ) ♦  2  fi £  sine  -2 12  n  cos 6  (24) 

ay(  =  ri  ♦(fl2xv  +  w*n20pcx)  sin8  ♦  (v  ♦  2ft£  -  fl2v 

ft2x  p)  cos  0  -  ft2  n  cos20  ♦  j  fl2  c  sin  2  0  (25) 
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.  I BSJWP 


=  c 


-(V  ♦  2  -  ft2v 


a* 


fl2  B^xJcose 


2  2  12- 
8  ;  sin  0  +  2  G  n  sin  20 


(26) 


where  (E.n.C)  is  a  general  point  relative  to  the  deformed  blade  coordinate  system.  It 
may  be  noted  that  the  gravitational  contribution  to  these  accelerations  has  been 

neglected. 


Flapping  Pendulum 

The  blade-mounted  pendulum  absorber  is  shown  in  figure  3.  The  pendulum  after 
deformation  is  shown  in  figure  4.  The  deformed  coordinates  of  the  pendulum  hinge  and 
mass  can  be  related  from  figure  4  and  are  given  by 


r  =  l  cos  A 
P 


n 


P 


nA 


r  =  r  .  ♦  l  sin  A 

*  p  A 


(27) 


where 

l  =  the  pendulum  arm  length 
A  =  rotation  about  (0,  nA,  CA)  in  the  x*  z'  -  plane 


|  =  0(6)  ;  A  =0(1) 


Differentiating  equations  (27)  with  respect  to  time  yields 


r  =  -l  A  sin  A  (28a) 

P 

r  s  -lAsinA-lA^  cos  A  (28b) 

P 


IS 


0 


(28c) 


n 

p 

=  0 

(28c) 

% 

=  0 

(  2 8d ) 

s 

a  4  A  cos  A 

(28e) 

•  2 

a  4(  A  cos  A  -  A  sin  A) 

( 2  8f  ) 

Substituting  the  above  expressions  into  equations  (24)  to  (26),  the  acceleration 
components  to  the  first  order  (second  and  higher  order  terms  are  discarded)  for  the 
flapping  pendulum  become 

a^,  =  -  Jl  A  sin  A  -4  A  2  cos  A  -211  vA  -  fl2xA 

xp 

-  fl2  4  cos  A  +  28  4  A  cos  A  sin  0A  (2  9) 

V  =  '“^A  VA*«A*  “V'*1  Sl"  ®A 

2  2 

+  (vA  -2  (1  4  Asin  A- fl  vA  +  fl  xA  yA)  cos  0  A 
2  2  12 

-  fl  nAcos  0  A  ♦  j  fl  (  c  A  ♦  t  sin  A)  sin  2  0A 

(30) 

•2  . 

a  ,  =  4  A  cos  A  -  4  A  sin  A 

P 

-(vA-2fl  4  AsinA-fl  vA+  0  xAwA>sin0A 

>(n2xAvA  +  wA+  fl2  6 pj.  XA)  cos  0  A 

-  fl2(cA*  4  sin  A)sin2eA*^fl2nA»in2  0A 

(31) 
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These  acceleration  components  will  be  used  to  express  the  inertial  reaction  force  and 
moment  at  the  pendulum  hinge. 

The  total  inertial  reaction  force  acting  on  the  pendulum  mass,  M,  can  be  wri$t$n  in 
terms  of  the  above  acceleration  components  whose  associated  unit  vectors  are  T',  j',k'  as 


I7  = 


-M  (a  ,  1*  + 
P 


a  .  T  ♦  a  ,  k') 
yP  P 


(32) 


The  position  vector  of  the  pendulum  mass  with  respect  to  the  hinge  point  is 


=  J,  (cos  A*'  +  sin  Ak') 


(33) 


and  the  inertial  moment  reaction  about  the  hinge  point  is 


M  =  rxF 


or 


-Ml 


cos  A  0 


sin  A 


ax’  ay>  az’ 

P  7  P  P 


M 


Mia,  sin  Ai' 
YP 


♦Ml  (a  ,  cosA-a  ,  sinA)j' 
P  P 


-M  l  a  ,  cos  A  k’ 
yP 


(34) 


In  addition  to  the  inertial  reaction  moment,  there  is  a  moment  due  to  damping  at 
the  pendulum  hinge.  This  damping  will  be  presumed  to  be  of  the  linear  viscous  type. 
The  corresponding  moment  reaction  about  the  hinge  can  be  represented  as 
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Ma  =  CNA? 


(35) 


where 


N  =  steady-state  normal  force  on  the  bearing  due  to  FAx,  and  F^, 

C  =  pendulum  hinge  friction  coefficient 

Since  the  flapping  pendulum  is  free  to  rotate  about  the  y'-axis,  the  y’  component  of 
the  inertial  moment  reaction  is  in  equilibrium  with  the  damping  moment.  Using 
equations  (34)  and  (35),  the  equilibrium  equation  can  be  written  as 

Mi,  (a,  sin  A -a,  cos  A)  -CN  A  =  0  (36) 

x  z 

p  p 


Equation  (36)  provides  the  pendulum  equation  of  motion.  The  other  two  components  of 
the  inertial  moment,  the  damping  moment,  and  all  three  components  of  the  force 
reaction  provide  the  pendulum  reactions  which  act  on  the  blade  at  the  hinge  point  as 


•Max'_ 


-May 


(37) 


“Ay- 


Mi  a  ,  sin  A 
XP 


=  -CN  A 


(40) 


(41) 


MAz'  =  -Mia  ,  cos  A 
yP 
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Substituting  for  the  acceleration  components  in  equation  (36),  the  pendulum 
equation  can  be  written  as 

....  2  2 

l  A  -  (vA-fl  VA  +  ft  xA  yA)  sin0A  cos  A 

2  ..  2 

+  (ft  XA  vA+wA  +  n  6p  XA*C0S  »aCOsA 

'  ■  “  *C 

2  2 

+  ft  8,  cos  0A  cos  A  sin  A 
2 

+  (2  ftvA  +  ft  xA)  sin  A 
2  2  12 

+  (-  ft  cA  sin  0  A  +  ^  ^  nA  sin  2  0^)  cos  A  =  -CN  A/M  11  (43) 

It  may  be  noted  that  in  the  above  equation  the  rotary  inertia  of  the  pendulum  is 
neglected.  This  implies  that  the  pendulum  consits  of  mass,  M,  which  has  very  small 
dimensions  (point  mass)  in  comparison  with  the  pendulum  arm. 

The  nonlinear  equation  of  motion  given  by  equation  (43)  is  linearized  for  small 
angular  oscillations,  fi,  about  the  steady-state  equilibrium  angle  of  the  pendulum,  6  , 
where  0 

A  r  6  ♦  6q  (44) 

Since  the  perturbation  angle,  6  ,  will  be  fairly  small  as 

52  «  1 

the  following  approximations  can  be  used. 

cos  a  =  cos  6Q  -  6  sin  6  Q  (45) 

sin  6  =  sin  6  ♦  6 cos  6  (4  6) 

00 

The  steady  state  equilibrium  angle,  6  ,  is  obtained  from  equation  (43)  by  setting 
the  acceleration  and  velocity  terms  to  zero°(assuming  that  the  elastic  deformations  are 
zero).  This  process  yields 

xAtan  6o^c°s2eAlsin  6o=  CA  sin2  eA  -  I  \  sin  2  0a  "  6pc  XA  cos  ®A  U7) 


In  the  above  equation  0A  is  the  steady-state  value  of  the  pitch  angle  at  xA.  It  may  be 
noted  that  the  steady  equilibrium  angle,  SQ,  is  independent  of  the  rotor  speed.  This 


result  is  due  to  neglecting  the  gravitational  force  in  equation  (36).  The  effect  of  the 
gravitational  force  on  the  steady  equilibrium  angle  and  the  uncoupled  pendulum 
frequency  is  very  small  (ref.  7). 

Substituting  equations  (45)  and  (46)  into  equation  (43)  and  subtracting  the  steady- 
state  relation,  equation  (47),  from  the  resulting  equation,  the  following  equation  of 
motion  for  the  perturbation  pendulum  angle  is  obtained. 

“  *sji  «  *  n2  S<ST6o  *  1  cos2  ®A  cos2  50>*  2«vA  Sin  S„ 

*  t(  VA  *  V  COS  ®A  '  (vA  -  5)2  VA  *  “^A  >'A)sin  »a]  cos  «o 


+  2  flvA  6cos  -  6  [(ft  2  xA  \>A  +  wA)  cos  0  ^ 
-(VA  *  VA  +  fl2  XA  PA)  Sin  6a]  Sin  6o  =  0 


In  equation  (48)  the  6th  and  7th  terms  (underlined  terms)  are  nonlinear.  Comparison  of 
the  7th  term  with  the  5th  term  illustrates  that  this  nonlinear  term  is  of  order  6  tan  6 
compared  to  unity  for  the  linear  term.  Since  both  6  and  <Sq  will  be  fairly  small  then 
the  7th  term  can  be  rationally  neglected.  Comparison  of  tne  nonlinear  6th  term  with 
the  linear  4th  term  indicates  that  they  are  of  the  same  order  of  magnitude  since 


Both  of  these  terms  will  be  fairly  small.  In  the  interest  of  linearization  and  the  desire 
to  retain  an  explicit  damping  term,  the  6th  term  will  be  dropped  and  the  4th  term 
retained. 

The  steady-state  normal  force  on  the  bearing  can  be  written  as 

N  =  [FAx'*FAz']I/2  (*»> 

where  F.  „  and  are  the  steady-state  values  of  FA  „  and  FA  „  respectively.  It 

may  be  noted  that  rv  is  generally  a  time-dependent  force  but  is  approximated  here  by 
the  steady-state  value.  Substituting  for  the  acceleration  components  (eqs.  (29)  and  (31)) 
in  equations  (37)  and  (39)  and  retaining  only  the  steady-state  terms,  equation  (49) 
becomes 


N  :  Mil  (xA  +  l  cos  6 Q) 


It  may  be  noted  that  high  order  terms  are  neglected.  Substituting  equation  (50)  into 
equation  (48),  the  linear  pendulum  equation  of  motion  becomes 


4  6  +  C  ft2  (y  ♦  cos  6q)  6  +  +4  cos2  0  A  cos2  sjs 

2 

+  2  ft  vA  sin  6q  +  (wA  +  ft  XA  v cos  9  A  cos  6 
-  (v A  -  ft2  vA  +  ft2  x A  ma)  sin  0A  cos  6q  =  0  (51) 


The  uncoupled  pendulum  equation  can  be  obtained  from  equation  (51)  by  dropping 
the  blade  deflection  terms.  The  resulting  equation  can  be  written  as 


6+2  to)  6  +  (i)  6  =  0 

p  p 


(52) 


where 


=  Q  (fcosV~  +  COs2  6o  cos2  6  a)1/2  <  5  3 ) 


5  = 


C  fi(xA/ 4  +cos6Q) 


2ll 


4  cos  6 


2  .  .  2  ]l/2 

+  cos  6q  sin  0,  r 


•> 


(54) 


It  may  be  noted  that  the  pendulum  hinge  coefficient,  C,  in  equation  (54;  can  not 
normally  be  specified,  but  representative  values  of  the  damping  ratio,  5  ,  can  be 
estimated. 

The  inertial  reaction  loads  applied  to  the  blade  due  to  the  pendulum  will  be 
considered  to  act  at  the  elastic  axis  pendulum  hinge  location,  (xA,0,0).  These  reactions 
will  be  treated  as  discontinuities  in  the  tension,  shear,  and  bending  moment  distributions 
along  the  blade.  The  values  of  these  discontinuities  can  be  written  in  terms  of  the 
component  reactions  as 
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! 


A  M  , 
x' 

=  MAx’  + 

nA  FAz*  "  Z  A  FAy' 

(55) 

4My, 

=  v+ 

C  A  FAx' 

(56) 

AMz, 

=  MAz-“ 

nA  FAx' 

(57) 

AVx> 

=  FAx' 

(58) 

AV 

=  FAy’ 

(59) 

AVZ, 

=  FAz' 

(60) 

Substitution  of  equations  (29)  to  (31)  and  equations  (37)  to  (42)  into  the  above  equations 
yields 

AMx,/M  =  ( £  sin  A  +  £*>  D  °2  XA  VA  +  *A  +  n  2  &pcXA*  sin  0A 

•  2  2 

+  (vA  -2  ft  l  Asin  A  -  ft  v  A  +  ft  xA  cos  0  A 

2  2  12  ~\ 

-11  nA  cos  0  A  +  2  n  (  CA  +  Sin  A)  sin  2  6  A  J 

-  tia  Qua  cos  A-  f,  A2  sin  A-(vA-2£l  f.  Asin  A 

2  2  2  " 

-  n  vA  +  n  xA  uA)  sin  0  A  +  ( £1  xA  vA  +  wA 

+  fi2BpcXA)cOS0A-n2  C  a  sir*2  9  A 

-  n2£sin  Asin2  ©A  +  j  0  2  nA  sin  2  0  A  ]  (61) 


/  • 
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•  2  2 

A  My,/M  =  cA  [i  A  sin  A  +  A  A  cos  A+  2ft  vA  ♦  ft  xA 


♦  ft2  A  cos  A  -2ft  £  Acos  Asin  ©A] 

-C  ft2  (xA  +  A  cos  A)  A 

2 

A  Mz,/M  =  -  nA  [*<A  sin  A+  i  A  cos  A+ 2ftvA 

+  ft2xA+  ft2  £cos  A-2ft  £  Acos  Asin  ©A] 

-  Acos  A  [(  ft2xA  vA  +  wA  +  ft2  B^x^  sin  0A 

+  (vA  -2  ft  A  A  sin  A-  ft2vA  +  ft2xA  pA)  cos  0A 

-ft2  nA  cos2  eA  ♦  j  ft  2(C  A  ♦  «•  sin  A)  sin  2  0A] 
..  .  2 

A  Vx,/M  =  £  A  sin  A  +  £  A  cos  A  +  2  ft  v  A 

+  ft  2  xA  +  ft2  £  cos  A-2  ft  £  Acos  Asin  0A 
A  Vy,/M  =  -<n2xA  vA  *  wA+  ft20pcxA)sin6A 

-  (v'A  -2  ft  £  Asin  A  -  ft2vA  +  ft2  xA  uA>  cos  0A 
♦  ft  2  nA  cos2  0  A  \  ft  2(  £  A  ♦  £  sin  A)  sin  2  0A 


(62) 


(63) 


(64) 


(65) 
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A  Vz^M  =  -4  Acos  A+  t  A  sin  A  +  (v‘A  -2  fl  i  A  sin  A 

-n2vA.  a2xA  MA)sineA-(a2xAvA«wA 

2  2  2 
+  fl  8  ^  xA)  c°s0  a  +  £1  zA  sin  0A 

+  SI  ^  l  sin  Asin2  0  A  ft  2  nA  sin  2  0  A  (66) 


where  M  is  the  pendulum  mass. 

Using  the  adopted  ordering  scheme  given  by  Table  1  and  equation  (28),  the  orders  of 
magnitude  of  the  above  expressions  are  as  follow 


AM  ,/M 
xw 

~ 

CXe2)  ' 

1 

\ 

A  MyyM 

- 

CKe)  ! 

>  (67) 

AMz,/M 

~ 

CKe)  j 

AVx,/M 

~ 

cxi)  y 

| 

AVyI/M 

~ 

CXe)  j 

i  (68) 

AVz,/M 

~ 

CXe)  ) 

1 

For  inclusion  in  the  blade  equations,  the  shear  and  moment  reactions  must  be 
transformed  into  the  undeformed  x,  y,  z  system  using  the  transformation  matrix  of 
Appendix  B.  The  transformations  consistent  with  the  above  orders  of  magnitude  are 

A  Mx/M  =  (  AMx,/M)  -  (A  My,/M)  (  pA  cos  0A  + 

vA  sin©  A 

) 

+  (  AMz,/M)  (  pA  sin  0A  -  vA  cos  0  A)  (69) 

A  My/M  =  (  AMy,/M)  (cos  0A  -  $A  sin  0  A)  -  ( A  Mz,/M)  (sin  0A  +  $A  cos  0A)  (7  0) 

AMz/M  =  (  AMy,/M)  (sin  0A  +  $A  cos  0A)  +  (  AMz,/M)  (cos  0a  -  $A  sin  ©A)  (71) 
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2 

where  only  O  ( t  )  terms  and  larger  are  retained.  Similarly 


AVx/M  = 

a  vxjm 

(72) 

AV  /M  = 

y 

(AV^M) 

yA  + 

( A  Vy,/M)  cos  0A  -(A  Vz,/M)  sin  0 A 

(73) 

AVz/M  = 

(avx;m) 

VA  + 

(  AVyI/M)  sin  0A  +  (A  Vz,/M)  cos  0A 

(74) 

where  only  O  ( e )  terms  and  larger  are  retained. 

Substituting  equations  (61)  to  (66)  into  equations  (69)  to  (74)  and  retaining  only  the 
linear  time  dependent  terms  (to  be  compatible  with  the  completed  time  dependent 
blade  equations),  the  shear  and  moment  reactions  in  the  undeformed  preconed  system 
are  obtained  after  making  use  of  equations  (44)  to  (46). 

A  Mx/M  =  [(4sin6Q  +  £A>  sin  0A  -  nAcos0A]wA 

+  [(4  sin  6o  +  tA)  cos  0A  +  nA  sin  0A]  [vA  -211  l  6  sin  «o 

2  2  2  ~\ 

-ft  vA+  ft  xA  pA  +  ft  4  6  cos  6q  sin  0A  J 

+  46  cos  6q  [  ft  2  8  ^  x  A  sin  0  A  -  ft  2  nA  cos2  9  A 

+  ft2(  5  a  +  4  sin  sin  ®acos  9A  3  "  nA4  ^  cos  5o 
2 

"  0  XA  ^  nA  Sin  ®A  +  ^A  C0S  0A*  mA 

+  4sin  6Qft2xA  vAsin  0A  (7  5) 
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1 


A  My/M 


A  Mz/M 


(nAsin  0A+  CA  cos  eA  >[*  6  sin  +  2ft  vA 

-  ft  2£  6  sin  6  -2  ft  l  6  cos  6^  sin  8,  1 
o  o  a J 

+  4cos6osin  8A  [( ft2  xA  vA  ♦  wA)  sin  0A 
+  (vA  -  2ft£6sin  6q  -  ft2  vA  *  ft2  xA  pA)  cos  0A 
♦  ft2£  6  cos  6q  sin  0A  cos  0A  ]  -ft  2  £  6  sin  6Q  sin  0A 

'  [VXASin9A-nACOs2  ®  A  +  (  ^ A 
+  £  sin  6q)  sin  0A  cos  8  ^  J  ♦  ft2  xA  (  nA  cos  0A 


-  C  A  sin  0A)  4>a  -C  ft  6(xa  +  £  cos  6q)  cos  0a 


(-  nA  cos  0A  +  cA  sin  0A)  [£  6  sin  6q 

♦  2  ftv.  -  ft2  £  6  sin  5  -2ft£  5  cos  6  sin  8,1 

A  O  O  A- 

-  4cos  6oCOS  eA  [(  XA  VA  +  WA>  sin  eA 

+  (vA  -2ft  £  6  sin  6q  -  ft2  vA  +  ft2  xA  VA)  cos  ©A 

+  ft 2£  5  cos  6q  sin  0A  cos  0AJ 
+  ft  2  £  6  sin  6Q  cos  0A  f  0  xA  sin  0A 

-  nA  cos2  0A  +  (  CA  +  sin  6q)  sin  8A  cos  0A  J 

2 

♦  ft  xA(  nA  sin  0A+  CA  cos  0  A)  4»a 


-C  ft  6(x,  ♦  £cos  5  )  sin  0, 

A  0  A 


(77) 
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A  Vx/M  =  1 6  sin  6Q  +  2  ft  vA  -  ft2£  6  sin  6Q  -  2  ft  £  ficos  6Q  sin  0A  (78) 


A  Vy/M  =  -vA  ♦  2  ft  i  6  sin  6Q  +  ft2  vA 


2 

-ft  £  6  cos  6q  sin  0A  +  £  6  cos  6q  sin  0  a  (79) 


A  Vz/M  =  - w  A  -  £  6  cos  6q  cos  0a  (80) 


Since  these  discontinuous  load  reactions  contain  the  time  dependent  pendulum 
angle,  it  must  be  related  to  the  elastic  displacements  for  its  elimination.  This  can  be 
accomplished  by  using  the  linearized  pendulum  equation  (eq.  (51) .  By  seeking  the 
simple  harmonic  motion  of  the  form 


6  (t) 


6eiut 


(81) 


and  corresponding  harmonic  solutions  for  the  state  variables,  the  pendulum  deflection 
can  be  written  as 


6  =  VA  6v  +  wA  6w+  “A  6y  +  Vv  (82) 

where 

6V  =  “  [(  ♦  <*>  cos  6Q  sin  0A  ♦  2  i  ft  a)  sin  6Q  J  /a  (83) 

2 

6W  =  w  cos  «Q  cos  0A/a  (84) 

5y  =  fl2xAcos  6osin  0A/a  (85) 

\  =  -  02  xA  cos  6q  cos  0A/a  (86) 
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a  =  A  (uip-w2)+  iCn2w  (-£  +  cos  6Q) 


(87) 


where  C  can  be  determined  from  the  estimated  damping  ratio  of  equation  (54)  as 


2  +  cos2  60  co$2  0A)] 

_ o  _ 

fl  (x^/A  +  cos  <SQ) 


The  discontinuous  load  reactions  associated  with  the  presence  of  the  pendulum  are 
used  to  relate  the  bending  moment  and  shear  force  state  variables  on  the  right  side 
(outboard)  of  x .  to  those  on  the  left  side  (inboard).  These  state  variables  are  indicated 
in  figure  5.  The  state  variables  are  related  as 


mr  , 

ml  -  am 

vR  = 

> 

< 

i 

-i 

> 

X 

X  X 

X 

X  X 

mr  = 

ml -  am 

VR  = 

> 

< 

1 

_l 

> 

y 

y  y 

y 

y  y 

mr  = 

ml-  am 

Z  Z 

vR  - 

z 

vL  -  AV 

Z  Z 

All  other  state  variables  must  be  continuous  across  the  point  x . . 

Substituting  for  the  shear  and  moment  reactions,  given  Dy  equations  (75)to  (80), 
with  the  explicit  time  dependency,  the  above  relations  can  be  described  by 

jzR)  =  [p]  jzLj  (89) 

where  £p3  is  a  "point"  transfer  matrix.  The  elements  of  [p[]  are  all  zero  except  for 


p..  = 

1 

(i  = 

j) 

P61  = 

m[ 

( A  sin  5^  + 

CA)  cos 

0A 

P62  = 

m| 

( 1  sin  6  + 

-  0 

CA)  sin 

0A 

P64  = 

-M 

i  sin  6  cos 
o 

9a  fl2 

XA  ' 

P65  = 

-M 

Asin  &  sin  0.  fi2  x 
O  n 

A  ' 

W  m  m 

v  x  6 


I 


-M  £(  nA  sin  0^  ♦  cA  cos  0A)  2iu>  ft 

2  2  1 

- £  cos  ^  sin  0A  cos  0A  (w  ft  )J  -  6y  M 

2  2 

M  a)  £  cos  ^  sin  0A-  Myfi 

2 

-M  (  nA  cos  eA  -  5  A  sin  0A)  ft  xA 
2 

-M  ft  x .  £  cos  6  sin  0,  cos  0.-6  M  , 
A  o  A  A  p  y{ 

-M  ft  2  x  »  £  cos  6  sin2  0.-6  M  r 
A  o  A  v  y  6 

-M  ^£  cos  6^  cos2  ©A  ( w2  +  ft2) 

-<nAcos  0A-  cAsin  8a)  2i(ua]  -  «yMz{ 

2 

-Ma>  l  cos  6o  cos  0A  sin  0A  *  Mz  g 

2 

-M  (  nA  sin  0A  +  cA  cos  0A)  ft  xA 

M  ft2  x .  £  cos  6  cos2  0.-6  M  . 

A  o  A  y  z6 

2 

M  ft  x »  £  cos  6  cos  0 .  sin  0.-6  M  . 

A  o  A  A  v  z  6 

-2  M  iuft  -  6  V  . 

v  x  6 


95  = 

-  6  V  , 
v  x  6 

10,1  = 

-M (w 2 ♦  n2).  «vvy6 

10,2  = 

-6  V  . 
w  y  6 

10,4  = 

-6  V  . 
y  y  6 

io,  y 

-  6  V  . 
v  y  6 

n,is 

•5  V  . 
v  z  6 

11,2s 

-Min2  -  6  V  . 

w  z  6 

ii,*s 

-  6  V  * 
y  z  o 

11,5  = 

-  6  V  . 
v  z  6 

where 


My  6  =  “M  |  ^  nA  sin  SA  *  ^A  COS  ®A^  1  sin  ^  +  ^  ^ 

+  £( £  sin  6q  +  Zp)  sin  0A  -  nA  cos  8^ 

2 

•fl  £  sin  6q  sin  &A  cos  8A  +  (x^  8  sin  &Q 
-  £ cos2  6q  cos  0 p)  H  2  £  sin2  &A  +  £ ( £  sin  6q 
+  CA)  cos  0A  +  nA  sin  0AJ  2ia)ft  £  cos  6q  sin  0A 
+  iC  Sl2u)  (x^  +  £cos  <SQ)  cos  8a  j 

!2  2 

(  nA  cos  0  a  -  cA  sin  eA)  £  sin  6q  <ui  +  fi  ) 

♦  [(£  sin  <5o  +  CA)  sin  0A  -  nA  cos  8A] 

2  2 

•  ft  £  sin  6  cos  6.  +(x4  S  sin  6 
o  A  A  pc  o 

-  £  cos2  5q  cos  0A>  cos  0A  sin  @A  +  nA  cos  8^ 
-  CA  sin  0A)  sin  0A  +  £  sin  <5q  cos2  @AJ 
•  2  i  wft  £  cos  6q  -iCfi  2u>  (xA  +  *•  cos  6  )  sin  ^A  J 


3U 


( 9  i ) 


(92) 


VX  s 

=  -M 

mm 

|_4  sin  (w  +  ft  )  ♦  2  i  u>  ft  4  cos  sin  8^  ] 

(93) 

vy  s 

=  -M 

f  4  cos  6  sin  0  .  (a)  ^  +  ft^)  -2  iu>  0  4  sin  6  1 

L-  o  A  o-» 

(94) 

Vz  6 

=  Mco 

2 

4  cos  6qcos  0a 

(95) 

Lead-Lag  Pendulum 

The  lead-lag  pendulum  after  deformation  is  shown  in  figure  6.  The  equations  of 
motion  for  the  lead-lag  pendulum  can  be  derived  by  following  the  procedure  used  in  the 
flapping  pendulum  case.  This  process  yields  the  following  equations  for  the  lead-lag 
pendulum: 

The  relation  between  the  pendulum  mass  and  hinge  in  deformed  coordinates 

£p  =  4  cos  A 

H  =  n.  +  4  sin  A 

P  A 

«p  =  «A 

Pendulum  angle 

A  =  6  +  6 
o 

where 

6q  =  steady-state  equilibrium  angle 
6  =  time  dependent  perturbation  angle 
The  equation  for  determination  of  the  steady-state  equilibrium  angle,  becomes 

2  1  2 
xA  tan  6q  ♦  4  sin  8^  sin  6q  =  -  ^  £ A  sin  2  0A  +  nA  cos  8^- ® pc  XA  sin  6 A 

The  pendulum  deflection  in  terms  of  the  blade  deflections 


5 


6  ♦  w. 
v  A 


6  + 
w 


VA 


6  + 
y 


VA 


6 

v 


(99) 
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where 


6 


6 


6 


6 


v 


w 


M 


v 


a 


to 


2 

P 


C 


6 


to 2)  cos  6q  cos  0a -2i  to  ft  sin  6q  ]/a 


=  0) 


2 

cos  6q  sin  0^/a 


2 

-ft  xA  cos  6q  cos  0A/a 


2 

-  ft  xA  cos  6q  sin  0A/a 


„  ,  2  2>  .  2  /  A  *  . 

H  C  ol>  -  to  if  iC!2  to  (—  +  cos  5  ) 
p  l  o 


=4 


XA  2  ,  .  2  „ 

- —  +  cos  6  sin  0 

cos  $Q  o 

x 


*) 


A  2^.2.  \l/2 

— r-  +  cos  6  sin  0  . 


2  - °~T~  +  cos*  sin*-  8  a  i 

^yjlcos  6q  o  Ay 

ft  (x./i.  +~  cos  5  ) 

A  O 


6eia,t 


The  nonzero  elements  of  the  point  transfer  matrix  become 


p..  = 

‘1 

1 

(i  =  j) 

P61  = 

M 

[cAc°s  eA  + 

(  nA  +  i  sin  6q)  -i.i 

i  0A  ]  (to  2  ♦  ft2)  - 

P62  = 

M 

fcAsin  eA-( 

nA  ♦  £  sin  6q)  cos 

0*  1  to  -  6  M  . 
A  J  w  x  6 

P64 

-M  *sin  ^  sin  0A  fl2  xA  -  6  ^  Mx  5 

(100) 

(101) 

(102) 

(103) 

(104) 

(105) 

(106) 

(107) 
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P65 

M  2,  sin  6^  cos  0 . 

0  A 

“2  *A 

-  6  Mv 

V  x  6 

11 

-M  [(  nA  sin  0A  + 

X,  A  cos  0A)  2iwft 

+  2,  cos  6 

2  2  ' 
•  sin  ©A(  to  +  ft  ) 

]-*, 

M  , 
y6 

P72  " 

M2,  cos  5  cos2  0 
0 

A“2- 

6M. 
w  y  0 

P73 

-M  (  nA  cos  ©A  -  i 

'-k  sin 

V  “2«a 

P74  = 

M  l  cos  SQ  cos  0A 

sin  9a 

“2xA'  5u 

M  x 

y  6 

P75  =  -M  2,  cos  6o  cos2  eA  n2xA  -  6V  My  6 
Pgj  =  -M  £2,  cos  6q  sin2  ©A( co 2  +  fl2) 

-(  nA  c°s  0A-  cAsin  eA)2ioifi]  -6vmz6 

Pg2  =  Mu,2&  cos  6osin  6Acos  0A-  6w  Mz  6 

2 

Pg3  =  -M(nAsin6A+  cAcos0A)JJ  xA 
Pg4=  Mn2xAicos6osin20A-6yMz6 

P85  =  'Mn2xA  1  cos  6o  Sin  6  A  cos  0A  '  6v  Mz  6 
P91  =  -2Miu,n-6vVx6 

P92  =  '  5w  Vx  6 

P94  =  "Vx* 

P95  -6v  Vx  6 


cos  0  . 
o  A 
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p  = 

10,1 

-M(u) 2  +  ft2)  -  6  V 
v  ; 

p 

*10,2 

P 

*10,4 

•  SuVy  6 

p 

10,5 

-  6  V  . 
v  y  6 

pn,r 

-  6  V  . 
v  z  6 

PU,2  = 

-Mai2  -6  V  „ 

w  z  6 

PHA= 

-  6  V 
p  z  6 

pn,r 

-  5  V  . 
v  z  6 

where 

Mx  6  =  ~M  j  [C  A  cos  ®A  *  (  nA  +  1  sin  6o)  sin  6 A  ] 

•  (ft  cos  6  cos  0.  +  io)  sin  5  )  2  ft  l 

O  A  O 

-  I  C  A  cos  6  Q  (  ft2  -  a) 2)  ♦  S  pc  ^xA  cos  6q  cos  0A  ]  (10  8) 


+  [~  5  A  sin  0A  +  (  nA  +  4sin  V  cos  0  A  ]  n2jt  sln  60sin  0ACOS  6  A 
♦  (xA  Bpc  sin  6o  "  4cos2  5o  sin  9a)  2|t  cos2  9  a 

'  [(  n  A  +  4  sin  V  Sin  0  A  +  cos  9 a]  2i  4005  50cos  9  A 
-  iCft2u(xA  +  a  cos  5q)  sin  0A  |  (109) 
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Mz  6=  M  |  (  nA  cos  0A  -  sin  0A)  £  sin  6Q  (to  2  +  ft 2) 

2  2 

+  £  C A  sin  0 A  -{  n A  +  £  sin  5q)  cos  0  A  J  ft  £  sin  6Q  sin*  0  A 


-(x.  0  sin  6  -£  cos*  6  sin  0.)  ft*  £  cos  8  »  sin  6  . 
A  pc  o  oA  A  A 

-  [(  nA  cos  0A  -  5A  sin  0A)  cos  0A  -  £sin  6q  sin2  0A  ] 

•2iwft£cos  6Q -iC  ft2w  (xA  +  8,  cos  6q)  cos  8Aj 

(110) 

x  6 

-M  ^£  sin  6q  (  u/2  +  ft2)  -2iuj  ft  £  cos  6qcos  0  A  ] 

(Ill) 

y<5 

M  ^£  cos  6q  cos  0A  (a)  2  +  ft2)  +  2iuft  £  sin 

(112) 

z6 

M  w  2  £  cos  6  sin  0. 

O  A 

(113) 

AERODYNAMIC  FORCES 

This  section  is  concerned  with  the  calculation  of  the  aerodynamic  loads  on  a 
helicopter  rotor  blade  in  forward  flight.  Due  to  the  complexity  of  the  flow  field 
involved  and  the  interaction  of  different  effects  in  the  flow,  a  completely  general 
solution  to  the  problem  has  not  been  obtained.  There  are  a  large  number  of 
approximate  methods  available.  These  methods  range  in  complexity  from  simple  blade 
element  representations  to  lifting  surface  models  with  freely  distorted  vortex  wakes, 
with  associated  ranges  in  computational  expense,  accuracy  and  detail  of  the  solution. 
Commonly  used  is  an  approach  based  on  a  combination  of  simple  momentum  theory  and 
the  blade  element  description.  This  method  is  adopted  for  the  present  study.  The 
method  takes  no  note  of  the  reverse  flow  and  stall. 

The  aerodynamic  forces  are  formulated  from  two-dimensional,  incompressible, 
quasi-steady,  strip  theory  in  which  only  the  velocity  components  perpendicular  to  the 
spanwise  axis  of  the  deformed  blade  are  assumed  to  influence  the  aerodynamic  loading. 
Account  is  taken  of  the  varying  freestream  velocity  associated  with  the  rotating  blade 
by  employing  Greenberg's  extension  (ref.  11)  of  Theodorsen's  unsteady  theory  for 
determining  the  aerodynamic  lift  and  pitching  moment  acting  on  the  blade.  The 
resulting  expressions  are  specialized  to  the  case  of  quasi-steady  flow  by  setting 
Theodorsen's  circulation  function  to  unity.  Simple  momentum  theory  is  used  to 
calculate  the  steady  flow  induced  by  the  rotor.  The  expressions  for  the  aerodynamic 
loading  are  based  on  the  analysis  reported  in  reference  12. 
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In  Theodorsen's  theory  a  two-dimensional  airfoil  is  assumed  to  be  pivoted  about  an 
axis  which  may  be  distinct,  in  general,  from  the  aerodynamic  center.  The  airfoil  is 
pitched  at  an  angle  e  (t)  to  the  f reestream  flowing  at  constant  velocity  V.  The  airfoil  is 
vertically  displaced  with  velocity  h  (t)  as  shown  in  figure  7.  Greenberg  (ref.  11)  has 
extended  Theodorsen’s  theory  for  a  varying  free-stream  velocity  V(t).  The  lift  and 
moment  acting  on  an  element  section  of  the  blade  may  be  expressed  in  terms  of  the 
circulatory  and  noncirculatory  components  as 
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cos  a  =  Uy/U 


u  . 

Substituting  equations  (115a),  (115b)  and  (116)  into  equations  (117)  yields 


Fy,  =  ipac[u2-|upc-^u|]  | 

f  (118) 

Fz'  *  5  P  [-UPUT  -  S  UT-  !  -  (S>2  =]  ) 

From  equations  (115c)  and  (11 5d)  the  pitching  moment  can  be  written  as 

Mx-(=V-5pac(t)2  [Ut'e  -Vt‘0  (119) 


x* ' 


where 

Up,  UT  =  relative  velocity  components  in  the  deformed  coordinate  system  (fig.  7) 
e  =  angular  velocity  of  the  blade  about  the  x'  -  axis 

V(t)  Zyj  Uy  + 

The  assumed  orders  of  magnitude  of  the  quantities  used  in  deriving  the  aerodynamic 
loads  are 

un 

J  .  OM 


>  =  OU2) 

a 


=  0(e) 


=  O  (e ) 


e 


O  (e ) 


The  aerodynamic  forces  in  the  undeformed  coordinate  system  are  obtained  by  using 
the  transformation  given  by  equation  (A-4).  To  second  order,  the  aerodynamic  loads  can 
be  written  as 


Lv(=Fy>  *  <y-Fz.<9  ♦♦> 

L„  <=<y  =  Fy.«**'*Fz, 


In  equations  (118)  and  (119)  the  blade  aerodynamic  loads  are  expressed  in  terms  of 
Uj,  Up,  and  e  .  In  order  to  use  the  expressions  for  blade  aerodynamic  loads  in  the 
blade  equations  of  motion,  Uy,  Up,  and  e  must  be  expressed  in  terms  of  the  blade 
bending  and  torsion  deflections  v,  w,  and  $  .  The  blade  velocity  is  easily  expressed  in 
the  x,  y,  z  coordinate  system.  The  deformed  blade  coordinate  transformation,  which 
relates  the  x,y,  z  and  x',  y',  z'  coordinate  systems,  is  given  by  equation  (A-4).  This  is 
used  to  express  the  Uy  and  Up  velocity  components  in  terms  of  v,  w,  and  $  . 

The  total  flow  velocity  seen  by  a  point  on  the  elastic  axis  can  be  written  from 
reference  1 2  as 

$  =  (pfftR  cos  ip  +  ft  RA  -u  +  v)T 

♦  (  pc  w  -v  -  ft(x  +  u)  -  pj  Q  R  sin 


where 


+  (  A  R X  -  p.  RR  8  cos  ip  -w  -  8B  v)  k  (121) 
i  pc  pc 


V  cos  a 
_ s 


and  the  aerodynamic  velocity  components  are  shown  in  figure  8. 

The  tangential  and  perpendicular  velocity  components,  Uy  and  Up,  are  obtained 
from  equation  (121)  using  £tJ  from  equation  (A-4).  To  secona  order  in  the  dependent 
variables  they  have  the  form 
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Uy  =  M|  fiR  v'  cos  i|>  +  ftx+  HjflRsin^  +  v  (122) 

Up  =  PjfiR  w*  cos\|»- pij  HR  v'8  cos\ti 

-(  yf  ft  R  sin4»  +  fix)  (0  +  $)-v0 


+  u.  £1R  0  cos  to  -flR  A  +  w  +  £10  v  (123) 
x  pc  pc 

It  may  be  noted  that  the  nonlinear  terms  are  neglected. 

The  angular  velocity  about  the  x'  axis,  e  ,  is  now  considered.  It  can  be  regarded  as 
being  composed  of  three  parts:  the  first  part  due  to  the  rigid-body  angular  velocity  of 
the  hub  in  space,  the  second  part  due  to  the  control  inputs,  and  the  third  due  to  the 
angular  velocity  associated  with  the  elastic  deformations.  The  total  sectional  pitching 
velocity  can  be  written  as  (ref.  12.) 


e  =  ft  (0pc+  W  +  0lcsin  ,l,'eiscos^+  ♦  (12  4) 


It  may  be  noted  that  the  total  geometric  pitch  angle,  0  ,  is  given  by 


0  =  0  pt  +  0  c  '  8  1c  C0S  *  ‘  0  ls  s^  ♦  (12  5) 


where  0  is  the  built-in  twist  angle  (pretwist),  0  is  the  collective  pitch  angle 
measured  «  the  tip,  and  0.  and  8 .  are  the  first  harmonic  cyclic  pitch  components. 
The  built-in  twist  is  a  linearlinear  function  of  x  and  may  be  written  as 


V  *  -V"-*’ 


(126) 


The  induced  veloctiy,  v.,  is  calculated  by  equating  the  integral  thrust  to  the  thrust 
from  momentum  theory  using  the  relations 


v. 

i 


CjftR 


(127) 
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c 


T 


1 

2  4 

u  p(TR* 


Lw  dxdij) 


(128) 


where  L  is  the  steady-state  value  of  L  . 

Substituting  equations  (118)  and  (11$  into  equations  (120)  and  using  expressions  for 
Up*  given  by  equations  (122)  to  (124),  the  aerodynamic  loads  to  second  order 

can  be  written  as 


Lv  =  5Pac[*^x+  yfORsin\ji)0  +  2  (  p{  RR0  cos 


R  X)J  to  +  ^  p  ac  £-(fl  x  +  PjftRsinvp) 


•(  y{  flR  Bpc  cos  i|>-«RX)J  $  +  |  pac  £-(flx 


+  Mj  n»l  sim|> )  0 


2(yfnR0pccos\|i-ftRX)J 


•  RR  cos  tyv  +  ^  p  ac  £-(ft  x  +  UjORsin^) 

Q 

•(y  .flR  3  cos\|»  -fl  R X )  0  — —  (fl  x 
i  pc  3 

+  PfflR  sin4»)2  +  (  pf  f)R  0^  cos  i|»  -fiRX  )2  1 
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j  p  ac  £-(ft  x  +  y^  ft  R  sin  \|>)  ft  6  ^ J  v 

+  ^  p  ac  £2  (ft  x  +  y^  ftR  sin\|» )  0  -(  y  ^  ftR  8^  cos  if; 

-ft  R  X  )J  v  +  ^  p  ac  £  -( ft  x  +  y  f  ft  R  sin  41)  w  -  |  w  J 

+  ±pac  [-(ftx+  yjftRsimjj)2 

+  4  R  cos  \JiJ  4>  +  i>  p  ac  (ft  x 

+  y  ^  ftRsimjj)^  4)  +  ^  pac  £  2  ( ft  x 

+  Pj  ftR  sin  40  8  -(y  jftRg  pC  cos  t -ft  RX  )J 

•(  y  j  ft  R  cost); )  y  +  ^  p  ac  £-(  ft  x 

+  y^  ft  R  sin40  y^  ftR  cos  i|»  +  ^  ft  (ft  x 

+  y  j  ft  R  sini|> )  +  ^  y ^  ft2  R  sin  4>  J  v 

+  ijP  ac  yf  ft  R  cos  \|»J  v 

+  ^  p  ac  £-(  yf  ftR  8^  co$4>  -  ft  RX  )  +  (  yf  ft  R  sin  4» 

♦  ft  x)  0^]  (  ftx  ♦  yj  ft  R  sin\|>) 

+  ^pac^|(n8pc+|0)^Rx  +  MfflR  sini|» ) 

♦  |  Mf  A2  R  cos  \|)0  +|  y{  ft2  R  0^  sin  4>J 


(130) 
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M  ^  ^  p  ac  (|)2  £-w  +  Wj  ft2  R  cos  4*  ♦ 

+  2(«x+  yf  fl  Rsin*)  ♦  +  (fl2x>  2  pf  C 

-  pf  ft  R  cos  i|i  \) ] 

-  ^  p  ac  (|)2  [(fix*  yfRR  sim|» ) 

•  (  R  8__  +  20)  +  p.  n2  R  8  cos  \|/ 

pc  I 

+  MfA2  R  6pc  sin  4<] 

The  above  expressions  can  be  rewritten  in  matrix  form  as 
dependent  terms) 


[BCx.t)]  |z(x,t)}  ♦  [c(x,  «>]{  Z  (x,  t)} 

♦  ^  [^Bc(x,  t)  J  n  {z  (x,  t)|  cos  n  i|> 

+  £  [Cc  (x,  t)  J  n  |  Z  (x,  t)  J  cos  n  ♦ 

+  V  [b  (x,  t)l  {  Z  (x,  t) }  sin  n  'P 
n  =  1 L  J 

*^^Cs(x,  t)]n  jz(x,  t)j  sinn  * 

*  j2fcc  «}  n  cos  "  *  *  jfc,  {Ps  <x,}nsin 


2 

!  R  sin  )  v 


(131) 

(retaining  only  time- 

+  Ip  (x,  t)j  |  z  (Xf  t)| 


♦  (132) 
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i 


where 


is  the  state 


vector,  defined  by  equation  (1),  and 


WM 


ooooo 


-M. 


0  0 


(133) 


vec^T.:"0  *h'  e*ast*c-^orma,*0,'-independent  aerodynamic  load 

For  n  =  lj 

pc6  =  jpac(|)2n  [-2 «x  els+  y{  nR(et+ 

PcI0=  -  j  Pac  OR  [-2  epc  yfRRX  -<0t+ 


+  4  6pcMf 

OR  9  is-AO  X9lc] 

li 

p  ac  fi  R 

L'ePc  utn*'  ,*  »c)  u,n 

is  "Hr  *  |  Rf  RR>  ®lc] 

Ps6  = 

|  pac  (|)2r 

[2n»8Ic*2a6?cM1R] 

o 

II 

-  j  P  ac  RR 

r  cdo 

L‘2T  +  <0  t+  ec)  Rx 

-X  fl  x  0  js  + 

ispcUf2BR8lc] 

Psll  = 

-  j  pac  R  R 

LT6pc  R  U,  +  X  (1  yf  R 

*  2(V  ec)  wfnx-£jr  OR)  els 

+  T«<|)  eic] 
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For  n  =  2; 


Rcg  =  2  P  ac  ^  ^  ^2  Pf^R  ®lc^ 


PclO  = 

-  j  p  ac  ft  R  £ 

1  Cdo  2  _  D  1  0  2  2  n  0 

2  a  MffiR  +  28pcMffiR 

+  ^  X  M ^  ftR 

6ls  +  iBpc  Mx  0lc] 

Pcll  = 

-  ^  p  ac  ft  R  | 

(e  t  +  ec>Hf  n  R  +  Q  Pf x  0  is 

"§  Mf  n  9  lc 

] 

Ps6  = 

\  P  ac  (|)2  ft 

(y  yf n  r  6  is} 

PslO  = 

pac  ftR 

["  2  (e  t  *  ec}  9  pc  yf  R  R  +  29  pc  Mf n  x  els 

-i  X  Mfft  R 

•lc] 

Psll  = 

-  j  p  ac  ft  R  £ 

2  9pcwf  R  '  2  c  vf  fi  9  Is  'ft  y  f  x  6lc] 

For  n  =  3; 

PclO  =  Pac  nR  6pc  Mj  ^  R  0ls) 

RCU  =  _2PacflR(^y£fJR  8jc) 

Psio=  -ipacflR^e^yJftR  elc) 

psn  =  -|pacnR(iyf fiR  y 
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The  first  seven  terms  in  equation  (132)  represent  the  aerodynamic  loads  due  to 
elastic  deformations.  The  coefficients  of  the  corresponding  matrices  are  given  in 
Appendix  C. 

The  aerodynamic  loads  as  represented  by  equation  (132)  are  to  be  included  in  the 
blade  equations  of  motion.  It  may  be  recalled  that  the  solution  to  these  equations  is 
treated  as  individual  contributions  associated  with  each  harmonic  of  the  excitation 
loading.  Such  a  solution  may  be  written  in  a  complex  form  as,  equation  (22), 


|zn(x,  t)|  =  Re(|fn(x)J  e1Wt)  (134) 

where 


a)  =  nfi  (n  =  1,  2,  3)  (135) 

In  the  aerodynamic  loads  expression  (eq.  (132))  the  fourth  to  the  seventh  terms 
contain  functions  of  the  azimuth  angle  (rotational  speed).  Due  to  the  existence  of  these 
terms,  the  above  solution  (eq.  (134)  can  not  be  applied.  Instead,  the  state  vector  should 
be  written  as 


{z(x,t)}=  Re^J  {zn(x)|einiitj  (136) 

In  this  case  the  resulting  equations  must  be  solved  simultaneously  for  all  harmonics. 
Even  if  only  three  harmonics  are  considered  for  the  state  vector  (eq.  (136))  obtaining 
the  corresponding  solution  is  very  complicated.  An  approximate  method  is  used  to 
determine  the  aerodynamic  forces  by  neglecting  the  fourth  to  the  seventh  terms  in 
equation  (132).  These  terms  are  of  lower  order  than  the  corresponding  terms  in  the 
blade  equations  of  motion  and  thus  may  be  neglected  according  to  the  adopted  ordering 
scheme.  Although  the  first  three  terms  of  equation  (132)  may  also  be  neglected 
according  to  the  ordering  scheme,  they  are  retained  here.  These  aeroelastic  terms  will 
be  the  only  representation  of  coupling  between  the  elastic  deformations  and  the 
aerodynamic  loading. 

With  the  fourth  to  the  seventh  terms  in  equation  (132)  neglected  the  solution  given 
by  equation  (134)  can  be  applied.  Substituting  equation  (134)  into  equation  (132)  yields 
the  following  complex  notation  for  the  nth  harmonic  of  the  aerodynamic  load. 

{Fa}  *  [\]  N  *{F}  U37, 


where 


amplitude  vector  of  the  aerodynamic  load 


amplitude  vector  of  the  elastic -deformation-independent  aerodynamic  force 
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and 


M  =  {Pc}-‘{Ps} 


fA  1  = 

|*B  (x,  go  ) 

+  i  co  C  (x,  co  ) 

-u>2  rD(x,«n 

L  aJ 

L 

L  J 

L  J 

From  equation  (139)  and  Appendix  C,  the  nonzero  elements  of  J  are 


Aa62=  fPac(f)2^2 


Aa63  =  I  Pac  ^2ia)  ^ 


Aa65  =  2  P  aC  V  Q  X 


Aaio  2  =  -^‘Pacw  [-fix(  8C+  0t)+|  p{  HR  0u-2SlRxJ 


a  10,3 


a  10,5 


all, l 


a  1 1,2 


^  p  ac  R2  x  R  X 


\  pac(|  r2x  Mf  r  6ic>  r2  M2R2Bpc) 


\  p  ac  {-si2  X  0^) 

^  ip  ac  a)  ^2R  x  (  0c  ♦  0  ^)  -  RR  0js -f  R  R  X  J 

1  c  2 

•  2  P  ac  (-i  u>  R  x  ♦  j-  go  ) 


Aall  3  =  ~1  Pac(-n2x2-i  p2  R2  R2  ♦  ^  iuRx) 


al  1,4 


5  »ac(-n2x „IRelc-ia2ufR28pc) 


(  138) 

(139) 
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vall,5 


=  -  ^  pac  (|  fi2  x) 


The  induced  velocity,  v.,  is  obtained  by  substituting  the  steady-state  value  of  L 
from  equation  (130)  into  equations  (127)  and  (128).  The  resulting  equation  is 


w 


'i  =  Pf  tanas-X  = 


an(6  X  :  Cn} 

ITU 


(140) 


where 


P*  = 


V  cosa  s 
ft  R 


(14  1) 


1  be 

an  =  48  TtRa 

Cn  =  (4  +  6  8C  -(1  +  3  y^)  0t  R  -6  pf  8ls 


(142) 


+  |(38 


pc 


4 


6lc> 


(143) 


Equation  (140)  is  an  algebraic  equation  in  X  . 


t 


TRANSFER  MATRIX  METHOD 
System  Equations 

The  blade  equations  of  motion  are  given  in  a  matrix  form  by  equation  (23)  for  the 
nth  harmonic  of  frequency  u)  as 


jZ(x)J’  =  [ A(x)]  {  Z  (x)j  +  jFa(x)j  (144) 

The  boundary  conditions  corresponding  to  a  hingeless  rotor  blade  are  idealized  in  this 
analysis  as  being  cantilevered  at  the  hub  (x  =  0)  and  free  at  the  tip  (x  =  R). 
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Hub:  {  zd  <0) }  =  l°l  (145) 

Tip:  [zf(R)J  =  j°|  (  1  4  6 ) 

where 

{  Zd  (x) J  T  =  Lv  w  ♦  y  vj  (14  7) 

jzt(x)JT  =  K  M,  m;  ?,  f(  fj  0*8) 


Equation  (144)  represents  a  set  of  11  linear  first-order  ordinary  nonhomogeneous 
differential  equations.  It  may  be  noted  that  is  a  complex  matrix  due  to  inclusion 
of  the  linear  damping  terms.  The  desired  solution  is  obtained  by  the  fomulation  of  a 
field  transfer  matrix  which  relates  the  state  variables  between  any  two  points  along  the 
blade. 

Transfer  Matrix 

The  concept  of  the  transfer  matrix  has  been  thoroughly  discussed  by  Pestel  and 
Leckie  in  reference  13.  The  backward  transfer  matrix  [T(x)3  is  defined  by 

(  Z(x)j  =  [T  (x)]  jz(O)}  (149) 


which  relates  the  state  vector  at  two  stations. 

It  is  often  possible  to  determine  the  transfer  matrix  by  using  simple  statics.  Such 
techniques  are  described  in  reference  13.  These  techniques  prove  satisfactory  for 
simple  lumped  parameter  or  lower  order  uniform  continuous  systems,  Dut  become 
cumbersome  when  used  for  nonuniform  and  higher  order  systems  because  of  the  required 
algebraic  manipulation.  Murthy  (ref.  9)  has  presented  a  systematic  approach  which 
eliminates  much  of  the  algebra  and  results  directly  in  differential  equations  for  the 
elements  of  the  transfer  matrix.  The  state  vector  usually  satisfies  the  differential 
equation 

^  jz  (x) j  =  [a(x)]  |z(x)  J  (150) 


Differentiating  equation  (149)  with  respect  to  x  gives 
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T,  |ZH  =  T*  [T«]  |zw> 


(151) 


From  equav^n  (149)  it  is  obvious  that 

Jz(0)j  =  [t(x)J_1  jz(x)j  (152) 

Substituting  this  into  equation  (151)  the  following  equation  can  be  obtained 

£|z<x)J  -  £  0<x>]  [t(x)]-1  {z<x)J  (15  3) 


Equating  equations  (150)  and  (153)  gives 

[[A(x)]  [T(x)]  [T(x)]"1  |Z(x)|  =  (°j  (154) 

Since  equation  (154)  must  be  satisfied  for  all  values  of  x  and  all  values  of  Z  it  follows 
that 

[aw]  =  A  [tw]  [tw]-1 

Then  post  multiplying  both  sides  by  (x)}  yields 

£  [t  (x)J  =  [A(x)]  [t  (x)  ]  (155) 


Therefore,  the  transfer  matrix  is  given  directly  by  the  solution  to  equation  (155).  By 
letting  x  go  to  zero  in  equation  (149)  the  initial  condition  of  the  transfer  matrix 
becomes  the  identity  matrix. 

[ T  (0)]  =  M  (156) 

If  equation  (155)  is  solved  as  a  coupled  set  of  first-order  differential  equations,  then 
equation  (156)  provides  a  sufficient  number  of  initial  conditions. 

The  transfer  matrix  of  the  blade  without  a  pendulum  can  be  obtained  by  integration 
of  equation  (155)  together  with  the  initial  conditions  of  equation  (156).  The  tension 
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coefficient,  T,  appearing  in  matrix  [ A]  is  given  by 


?  R 

T  =  ft  1  J*  mxj  dxj  (157) 


When  a  pendulum  is  attached  to  the  blade  the  transfer  matrix  at  any  spanwise 
location  is  obtained  as  follows.  Let  x^  be  the  spanwise  location  of  the  pendulum  hinge: 

Case  1.  0  <  x  <  x  A 

The  transfer  matrix  in  this  case  is  obtained  by  the  direct  integration  of  equation 
(155)  together  with  the  initial  conditions  given  by  equation  (156)  similar  to  the  case  of 
the  blade  without  the  pendulum,  but  the  tension  coefficient  is  given  by  the  following 
equation  instead  of  the  one  given  by  equation  (157). 


9  /•  ? 

T  =  ft  J  mxjdXj+Mft  (x A  +  £  cos  6q)  ( 1  5 S ) 


Case  2.  x^  <  x  <  R 

Let  ^T.  ]  be  the  transfer  matrix  of  the  blade  at  x  =  x^;  this  matrix  can  be  obtained 
from  Case  l  as  described  above.  Recall  that  [  P]  (see  eq.  (89))  is  the  point  ti^nsfer 
matrix  of  the  pendulum.  Let  [T2  (x)]  be  the  transfer  matrix  of  the  blade  from  x^  to  x. 
This  matrix  is  obtained  by  the  integration  of  equation  (155)  with  the  initial  conditions 
given  by  equation  (156)  from  x to  x.  While  integrating  this  equation  the  tension 
coefficient  given  by  equation  (137)  should  be  used  instead  of  the  one  given  by  equation 
(158).  By  the  product  rule  of  transfer  matrices,  the  transfer  matrix  of  the  system  at 
any  spanwise  location  for  this  case  is  given  by 

[T<x>]  =  [t2<x>]  [p]  [t,]  (15  9) 


Once  the  system  transfer  matrix  is  obtained,  the  state  vector  at  any  spanwise  location 
can  be  determined. 

System  Dynamic  Response 

It  has  been  shown  how  the  transfer  matrix  can  be  used  to  solve  the  homogeneous 
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differential  equations.  An  application  of  such  a  solution  is  to  find  the  natural 
frequencies  and  mode  shapes  of  an  elastic  system  (ref.  14).  Once  they  are  known  it  is 
possible  to  solve  the  most  general  cases  of  forced  undamped  vibrations,  whether  they 
are  transient  or  steady-state  in  character,  either  for  discrete  or  continuous  systems, 
through  the  normal  mode  approach.  On  the  other  hand  the  steady-state  response  caused 
by  a  harmonic  excitation  is  more  readily  solved  with  the  aid  of  a  particular  integral  of 
the  nonhomogeneous  differential  equation  without  making  use  of  the  normal  modes  and 
natural  frequencies  (ref.  (13)). 

The  procedure  adopted  here  for  the  steady-state  response  of  forced  vibrations  with 
harmonic  excitation  is  usually  called  the  "extended"  transfer  matrix  approach  for 
undamped  systems  and  the  "complex  extended"  transfer  matrix  approach  for  damped 
systems.  The  procedure  is  valid  for  any  linear  damping  since  normal  modes  are  not 
introduced.  The  restrictions  are  that  the  response  is  steady  and  excitations  are 
harmonic.  With  the  knowledge  of  the  system's  response  to  harmonic  inputs,  the 
response  to  any  arbitrary  input  can  be  computed  by  synthesizing  the  arbitrary  forcing 
function  from  an  aggregate  of  infinitesimal  harmonic  forces. 

The  blade-pendulum  equation  of  motion  was  given  by  equation  (23)  as 

(z(x)J'  *  [a(x)]  jf(x)j  ♦  ( Fa  (x)  j  (160) 

By  definition  of  the  transfer  matrix,  the  solution  for  the  homogeneous  part  of  equation 
(160)  can  be  written  as 

KM  =  [t(x)]  l^(o)i  <i6i> 

The  particular  solution  of  equation  (160)  can  be  written  as 

[  Zp  <x>j  =  f  £t  (x,  s)J  jpa^j  ds  (162) 

where  [T  (x,s)3  is  the  transfer  matrix  of  the  system  from  s  to  x  (in  contrast  to  ^T  (x)], 
the  transfer  matrix  from  0  to  x).  By  the  product  rule  of  transfer  matrices  the  following 
relation  can  be  written 

[ T  (x)J  =  [t  (x,  s)]  [t(s)]  (163) 

From  equations  (162)  and  (163)  the  particular  solution  can  be  written  as 

J  Zp  (x)  J  =  [T  (x)]  /  [T  (s)]  -1  j  F(s)J  ds  (16  4) 

The  complete  solution  of  equation  (160)  can  be  obtained  by  adding  equations  (161)  and 
(164). 


{ Z  (x)  |  =  [  T  (x)  J  |  Z  (0)  j  +  [t  (x)]  /[T(s)]*1  {F(s)j  ds  (165) 


Equation  (165)  gives  the  amplitude  vector  of  the  simple  harmonic  response.  The  vector 
j  Z  (0)  |  is  not  completely  known,  but  it  can  be  determined  from  the  boundary 
conditions  of  the  system.  _ 

To  determine  the  vector  {  Z  (0)  J ,  first  rewrite  equation  (165)  as 


( zd (x)  i 

i 

'V*>  1  Vx>' 

( z‘f w ! 

P 

_TC(X)  ;  Td(x) 

/  - 


zd(0) 


Zf  (0) 


where 


■Ta  (X)  (  Tb  (xf 
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c 

■  sa  (S)  , 

Sb  (s)- 

(  Fd  (s)  ) 

+ 

- 1 - 

/ 

- + - 

\ - >  ds 

Tc  (x)  1  Td  (x)_ 

ISc(s)  ' 

Sd(s)_ 

(  Ff  (s)  ) 

(166) 


T 

= 

[v  w  <}»  \1  vj 

lz'.l T 

= 

I  M  M  M  V  V  V  I 

Lx  y  z  x  y  zj 

lFdl 

= 

H 

(167) 

IP  T 

= 

l-M.  0  0  0  -L  -L  1 

L  4>  v  wJ 

(  168) 

[s  (s)] 

= 

M" 

(169) 

The  boundary  conditions  corresponding  to  a 

hingeless  rotor  blade  are 

|v»l 

= 

PI 

( 170) 

1  \  (R>  | 

= 

PI 

(171) 

Substituting  equations  (167)  and  (170)  into  equation  (166)  yields 
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From  equation  (172) 
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(172) 


Substituting  the  boundary  condition  of  equation  (171)  into  equation  (173)  yields 

{ v°>}  =  -  [vr)]'‘  [Tc(R) :  Tc  HZ  Ip  Ml* 


(173) 


( 174) 


Equation  (172)  gives  the  response  state  vector  of  a  hingeless  rotor  blade  due  to  a 
simple  harmonic  load.  The  hub  constraint  reactions,  {  ZAO) }  ,  which  appear  in  equation 
(172)  can  be  obtained  from  equation  (174).  It  may  be  noted  that  the  transfer  matrix  and 
all  state  variables  are  complex  due  to  inclusion  of  the  linear  damping  terms  in  the 
system  equations  of  motion. 

The  transfer  matrix,  by  definition,  is  independent  of  the  boundary  conditions. 
Having  determined  the  transfer  matrix,  the  boundary  conditions  can  be  simply  imposed 
to  determine  the  response  state  vector  as  shown  above  for  the  case  of  a  hingeless  blade. 
So  all  the  problems  corresponding  to  various  boundary  conditions  can  be  analyzed  easily. 
Another  advantage  of  the  method  is  that  it  provides  for  the  continuous  integration  of 
the  equations.  In  the  present  study  the  Runge-Kutta  procedure  is  utilized  to  solve  the 
differential  equations  governing  the  elements  of  the  transfer  matrix  (eq.  (155)).  One 
distinguishing  advantage  of  the  transfer  matrix  is  that  the  hub  constraint  reactions  are 
obtained  as  a  direct  result  (eq.  (174)).  The  extended  transfer  matrix  approach  is  valid 
for  any  linear  damping  since  normal  modes  are  not  considered.  The  number  of  stations 
at  which  the  transfer  matrix  is  determined  depends  only  on  the  step  size  of  the 
integration.  The  incremental  step  size  can  be  selected  arbitrarily  and  then  the  result 
compared  with  that  obtained  for  a  smaller  increment.  In  this  way  an  optimal  value  for 
the  step  size  is  obtained. 


RESULTS  AND  DISCUSSION 


The  object  of  this  study  is  to  establish  a  comprehensive  analytical  design  procedure 
for  the  installation  of  simple  pendulums  on  the  blades  of  a  helicopter  rotor  to  suppress 
the  force  and  moment  reactions  at  the  hub  during  forward  flight  conditions.  This 
procedure  will  correspond  to  a  process  of  optimization  with  the  ultimate  object  of 
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minimizing  the  root  shear  and  moment.  The  minimization  of  these  reactions  will  be 
accomplished  by  a  systematic  variation  of  the  pendulum  parameters  (both  geometrical 
and  inertial).  In  this  section  numerical  results  are  presented  together  with  a  discussion 
of  their  significance.  The  numerical  results  presented  pertain  to  both  uniform  and 
nonuniform  hingeless  rotor  blades.  The  properties  of  these  blades  are  given  in  Appendix 
D. 

Simple  Flapping  Pendulum 

The  first  part  of  this  section  compares  the  root  reactions  from  this  study  with  the 
results  of  Murthy  (ref.  1 5)  for  an  undamped  analysis.  This  comparison  is  made  solely  to 
verify  the  correctness  of  the  computer  program.  The  remainder  of  this  section 
contains:  1)  a  detailed  discussion  of  the  numerical  results  for  a  4/Rev  concentrated 
load  at  the  tip;  and,  2)  the  results  for  excitation  by  a  harmonic  variation  of  spanwise 
airload  in  forward  flight. 

Analysis  verification.  -  By  dropping  the  damping  terms  in  the  blade  and  pendulum 
equations  of  motion  the  case  of  undamped  dynamic  response  can  be  recovered.  The 
undamped  root  reactions  of  a  uniform  rotor  configuration  (Appendix  D)  are  presented  in 
Table  2  together  with  the  results  of  reference  15.  The  reactions  with  a  properly  tuned 
simple  pendulum  are  also  shown  in  Table  2.  The  excitation  force  for  all  these  results  is 
a  concentrated  out-of-plane  simple  harmonic  load  at  the  tip  of  the  blade  with  a 
frequency  equal  to  4/Rev.  It  is  apparent  from  the  data  presented  in  Table  2  that  the 
present  analysis  provides  results  which  agree  with  the  work  of  reference  15  for  a 
uniform  rotor  blade  without  damping.  The  computer  code  has  also  been  verified  by 
duplicating  published  data  of  normal  vibration  modes  for  hingeless  and  articulated  rotor 
blades  (ref.  16). 

Concentrated  load  excitation.  -  The  force  vector  for  a  concentrated  out-of-plane 
load  can  be  written  as  (see  eq.  (168)) 


{f,}T=  Looooo  -fJ  (175) 

where  F  is  the  amplitude  of  the  applied  load.  The  concentrated  harmonic  force  at  the 
tip  has  tfie  following  characteristics: 

Magnitude . 2224  N 

Frequency . 4/REV 

Both  the  uniform  and  nonuniform  blades  are  considered  for  this  particular  excitation. 

The  root  reactions  of  the  uniform  blade  with  a  flapping  simple  pendulum  are 
computed.  It  may  be  noted  that  all  state  variables  are  complex  due  to  inclusion  of  the 
linear  damping  terms.  The  results  for  the  vertical  shear  amplitude  at  the  hub  are 
plotted  in  figure  9.  The  pendulum  weight  is  66.726  N  (10%  of  the  blade  total  weight). 
Damping  at  the  pendulum  hinge  is  neglected.  The  effect  of  this  damping  on  pendulum 
tuning  is  discussed  later  in  this  section. 

The  tuning  of  the  pendulum  is  changed  by  varying  the  spanwise  location  of  the 
pendulum  for  a  fixed  value  of  the  uncoupled  pendulum  frequency.  It  may  be  noted  from 
equation  (53)  that  the  uncoupled  natural  frequency  of  the  pendulum  is  mainly  dependent 
on  (1)  rotational  speed  of  the  blade,  (2)  spanwise  location  of  the  pendulum,  and  (3)  its 


TABLE  2.  -  CORRELATION  OF  UNDAMPED 
ROOT  REACTIONS 


Source 

M  , 
x’ 

Nm 

M  , 

y’ 

Nm 

Mz, 

Nm 

V 

N 

N 

No  Pendulum  Attached 

Reference  15 

-139 

1536 

-648 

-1054 

3256 

Present 

-146 

1535 

-648 

-1051 

3254 

With  Pendulum  Attached3 

Reference  15 

-5 

58 

-576 

-983 

-0.5 

Present 

-23 

53 

-606 

-973 

-2.8 

a  a>p  =  140  rad/s;  xA  =  1.651  m  ;  &A  =  0.2032  m 
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length.  The  simple  pendulum  considered  here  swings  in  the  flapping  plane  and  hence  can 
generate  significant  out-of-plane  shear  force  and  moment  at  its  point  of  connection 
with  the  blade.  These  facts  can  be  observed  from  the  results  of  figure  9.  The  root 
normalized  vertical  shear  (  the  shear  reaction  with  the  pendulum  attached  to  the  blade 
divided  by  that  value  without  the  pendulum)  is  altered  significantly  by  the  pendulum. 
For  instance,  when  the  pendulum  is  located  at  15.2%  of  the  span,  the  root  normalized 
vertical  shear  is  attenuated  to  2.9%  for  a  pendulum  frequency  equal  to  146  rad/sec 
(3.87/Rev),  while  it  is  amplified  to  2579%  for  a  pendulum  frequency  equal  to  160 
rad/sec  (4.24/Rev).  The  latter  case  shows  that  an  improperly  tuned  pendulum  might 
amplify  the  root  reactions  rather  than  attenuate  them,  and  the  system  would  be  better 
off  without  the  pendulum. 

From  the  results  of  figure  9,  it  can  be  observed  that  the  pendulum  of  frequency  146 
rad/sec  (3.87/Rev)  located  at  18%  of  the  span  is  a  proper  design  because  it  attenuates 
the  hub  vertical  shear  to  1.3%.  The  out-of-plane  moment  at  the  hub  is  reduced  to 
14.2%.  This  choice  of  pendulum  also  attenuates  the  in-plane  force  and  moment  at  the 
hub  to  4.2%  and  22.8%,  respectively. 

The  effect  of  pendulum  mass  on  the  vertical  reaction  is  illustrated  in  figure  10  as  a 
function  of  spanwise  location.  The  results  are  for  a  simple  flapping  pendulum  of 
frequency  146  rad/sec.  It  can  be  observed  from  these  results  that  if  the  weight  of  the 
pendulum  is  66.726  N  (at  x.  =  18%  span)  then  the  normalized  vertical  shear  at  the  hub 
is  1.3%.  If  the  weight  of  tne  pendulum  is  reduced  by  two  thirds  (located  at  18%  span) 
then  the  normalized  vertical  shear  at  the  hub  is  39.3%.  This  variation  of  the  pendulum 
weight  will  be  referred  to  as  the  weight  ratio,  where  in  this  instance  the  weight  ratio  is 
33%.  This  result  indicates  that  even  though  the  uncoupled  natural  frequency  of  the 
pendulum  is  independent  of  the  mass  of  the  pendulum,  a  proper  choice  of  the  mass  is 
also  required  to  generate  an  optimum  force  to  attenuate  the  root  reactions. 

It  can  be  observed  from  figure  9  that  the  minimum  normalized  vertical  shear  for 
any  pendulum  frequency  occurs  when  the  pendulum  is  located  near  18%  of  the  span. 
These  values  are  replotted  for  different  uncoupled  pendulum  frequencies  in  figure  11. 
Comparison  of  the  hub  vertical  shear  with  pendulum  absorber  to  the  hub  shear  with  no 
absorber  shows  the  pendulum  absorber  to  be  effective  in  reducing  hub  out-of-plane 
reaction  for  tunings  between  3.82  and  4.00/Rev.  It  can  be  recalled  that  the  optimum 
tuning  is  3.87/Rev.  The  steep  gradient  of  the  hub  shear  with  pendulum  tunings  beyond 
optimum  indicates  decreasing  effectiveness  of  the  pendulum  absorber  and  for  tunings 
above  4/Rev  the  pendulum  absorber  becomes  detrimental.  Figure  11  demonstrates  that 
in  the  case  of  an  off-tuned  pendulum,  it  is  apparent  that  the  result  can  be  either  a 
slight-to-moderate  degradation  in  pendulum  absorber  effectiveness  or  a  drastic  increase 
in  hub  reaction  depending  upon  the  direction  of  tuning  variation  from  the  optimum 
pendulum  natural  frequency. 

Similar  results  have  been  obtained  in  reference  4.  In  this  reference  a  semi- 
empirical  analysis,  using  in-flight  measured  flap  bending  moment  data,  was  developed  to 
study  the  pendulum  effectiveness  in  reducing  4/Rev  blade  root  vertical  shear.  The 
Boeing  Vertol  347  aircraft  was  used  for  the  flight  test.  It  is  a  tandem-rotor  helicopter 
with  fully  articulated  4  bladed  rotors.  The  vertical  4/Rev  pendulum  absorbers  were 
installed  on  all  blades  of  both  rotors  at  16%  blade  radius.  Comparison  of  the  predicted 
hub  motion  with  the  pendulum  absorber  showed  that  the  pendulum  absorber  is  effective 
for  tunings  between  3.4  and  4.0/Rev.  Test  results  indicated  that  a  3.95/Rev  is  an 
optimum  tuning.  For  this  optimum  tuning  the  vertical  hub  shears  were  reduced  by  80  to 
90%  on  both  forward  and  aft  rotors.  However,  the  vibratory  blade  root  chordwise 
bending  was  increased  by  25%  with  the  pendulum  installed. 

The  elastic  deflections  and  structural  load  reactions  for  the  present  uniform  blade 
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corresponding  to  the  optimum  flapping  pendulum  (3.87/Rev  at  18%)  are  shown  in  figure 
12  together  with  values  obtained  without  the  pendulum.  The  chordwise  and  flapwise 
deflections  are  normalized  with  respect  to  the  chord  (c=0.7m).  It  may  be  recalled  that 
all  the  deflections  and  forces  are  complex.  The  phase  angles  for  the  deflections  are 
found  to  be  either  0°  or  180°  as  presented  in  figure  12.  However,  the  phase  angles  for 
the  moments  and  forces  are  not  exactly  0°  or  180°  at  a  few  blade  stations,  and 
therefore  it  is  appropriate  to  plot  only  the  magnitude  of  these  forces  and  moments. 

It  can  be  observed  from  figure  12  that  the  pendulum  causes  a  spanwise 
redistribution  of  the  structural  loads  such  that  the  hub  reactions  are  attenuated.  Figure 
12  shows  that  the  optimum  pendulum  decreases  the  loads  inboard  of  the  pendulum  and 
increases  the  loads  outboard  of  the  pendulum.  It  can  also  be  observed  that  the  maxima 
of  the  deflections,  twisting  moment,  and  vertical  shear  are  higher  than  the 
corresponding  quantities  without  the  pendulum.  This  is  the  penalty  one  must  pay  for 
achieving  the  drastic  reductions  in  the  root  reactions.  The  vertical  shear  is 
characterized  by  a  large  change  at  the  pendulum  attachment  point  due  to  the  vertical 
inertial  reaction  of  the  flapping  pendulum  mass.  This  vertical  reaction,  which  is 
imposed  at  the  pendulum  attachment  point,  is  the  basis  of  the  blade-mounted  flapping 
pendulum  concept  which  has  the  purpose  of  attenuating  the  root  force  and 
moment  reactions. 

In  general,  a  properly  tuned  pendulum  can  attenuate  the  vibratory  loads  in  two 
ways:  (1)  by  eliminating  the  resonant  responses  of  the  blade  by  displacing  the  initial 
natural  frequencies  that  are  in  the  vicinity  of  the  excitation  frequency,  and  (2)  by 
generating  appropriate  forces  at  its  attachment  point  with  the  blade.  These  forces 
redistribute  the  loads  on  the  blade  so  that  the  root  forces  are  attenuated  as  discussed 
above. 


Figure  13  shows  the  hub  vertical  shear  for  pendulum  tuning  of  the  nonuniform 
blade.  The  blade  rotor  speed  and  precone  angle  are  equal  to  289  rpm  and  0.1  rad, 
respectively.  The  pendulum  weight  is  133.452  N  (13.8%  of  the  blade  total  weight).  The 
results  indicate  that  a  pendulum  with  uncoupled  natural  frequency  of  3.85/Rev  located 
at  22%  of  the  span  is  a  proper  design  because  it  attenuates  the  hub  vertical  shear  to 
2.2%  and  the  out-of-plane  moment  to  3.4%.  This  choice  of  pendulum  attenuates  the 
root  in-plane  force  to  47.6%.  However  the  root  in-plane  moment  is  amplified  by  75%. 
This  adverse  effect  can  be  reduced  by  offsetting  the  pendulum  hinge  (  n  ^ )  from  the 
elastic  axis.  The  effect  of  the  pendulum  hinge  offset  is  discussed  later  in  this  section. 

For  the  optimum  pendulum  (u  =  3.85/Rev  and  x^  =  22%  span)  the  root  forces  are 
computed  for  different  values  of  tne  pendulum  mass.  The  results  are  shown  in  figure 
14.  It  can  be  observed  that  a  50%  weight  ratio  (66.72  N)  is  an  optimum  choice.  For  this 
weight  the  normalized  vertical  shear  at  the  hub  is  practically  zero  as  compared  to  2.2% 
with  a  pendulum  of  weight  133.452  N.  This  demonstrates,  once  again,  that  a  proper 
choice  of  pendulum  mass  is  required  to  minimize  the  root  forces. 

In  the  above  analysis  the  damping  at  the  pendulum  hinge  is  neglected.  However,  it 
is  necessary  to  study  the  effect  of  the  pendulum  damping  on  the  root  reactions  once  the 
optimum  pendulum  tuning  is  established.  The  hub  shears  and  moments  are  computed  for 
different  damping  ratios,  and  these  results  are  presented  in  figure  15.  It  can  be 
observed  from  these  results  that  damping  of  the  pendulum  due  to  friction  in  the  pivot 
bearing  has  a  very  small  effect  on  the  hub  reactions. 

Once  the  optimum  tuning  configuration  and  mass  effectiveness  are  established,  a 
parametric  study  is  conducted.  The  parameters  to  be  varied  include  prepitch,  pretwist, 
precone  and  pendulum  hinge  offsets.  The  parametric  survey  is  conducted  for  the 
nonuniform  blade  with  the  optimum  tuned  pendulum  (  uj  =  3.85/Rev  and  x*  =  22% 
span).  The  pendulum  mass  is  133.452  N  for  this  study.  " 
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Figure  16  shows  the  effect  of  root  collective  Ditch  (prepitch)  on  the  hub  reactions. 
The  pretwist  angle  for  these  results  is  equal  to  -9  .  The  prepitch  has  a  small  effect  on 
the  hub  vertical  shear.  It  significantly  affects  the  hub  in-plane  force  and  moment, 
which  are  both  increased  by  increasing  the  root  collective  pitch.  The  out-of-plane 
moment  is  affected  moderately  by  changing  the  collective  pitch.  The  effect  of 
collective  pitch  on  the  twisting  moment  (not  shown  in  figure  16)  is  negligible. 

The  effect  of  the  built-in  twist  angle  (pretwist)  on  the  hub  reactions  of  the  blade- 
mounted  flapping  pendulum  is  shown  in  figure  17.  The  root  collective  pitch  for  these 
results  is  equal  to  13°.  Figure  17  shows  that  the  pretwist  has  a  small  effect  on  the  hub 
vertical  shear.  Unlike  the  effect  of  the  collective  pitch,  the  pretwist  has  a  significant 
effect  on  the  root  out-of-plane  moment.  However,  for  typical  rotor  blades  (negative 
pretwist)  the  root  out-of-plane  moment  is  still  attenuated  to  less  than  10%.  It  can  be 
observed  from  figure  17  that  the  in-plane  force  and  moment  are  both  increased  by 
decreasing  the  blade  pretwist.  The  pretwist  has  negligible  effect  on  the  twisting 
moment.  It  can  be  concluded  that  the  flapping  pendulum  effectiveness  in  attenuating 
the  root  out-of-plane  moment,  and  the  in-plane  force  and  moment  is  reduced  by 
decreasing  the  blade  pretwist. 

The  root  reactions  of  the  blade  are  computed  for  various  precone  angles,  and  these 
results  are  presented  in  figure  18.  It  can  be  observed  that  the  blade  precone  has  no 
significant  influence  on  the  hub  reactions.  Thus,  the  flapping  pendulum  effectiveness  in 
attenuating  the  blade  root  reactions  is  not  affected  by  changing  the  precone  angle. 

The  results  from  varying  the  pendulum  hinge  chordwise  offset,  n  are  shown  in 
figure  19.  It  may  be  recalled  that  the  optimum  flapping  pendulum  is  tuned  with  its 
hinge  located  at  the  shear  center  of  the  blade  cross  section  (  n*  =  0,  CA  =  0).  It  can  be 
argued  that  this  is  an  impracticai  location  since  the  shear  center  is  usually  inside  the 
airfoil  envelope.  The  pendulum  chordwise  hinge  offset  has  a  small  effect  on  both  the 
root  out-of -plane  force  and  moment  as  shown  in  figure  19.  Both  the  in-plane  force  and 
moment  are  decreased  by  increasing  the  chordwise  offset.  The  effect  of  chordwise 
offset  on  twisting  moment  at  the  hub  is  negligible. 

Figure  20  shows  the  root  forces  for  different  values  of  the  pendulum  hinge  vertical 
offset,  It  can  be  concluded  from  figure  20  that  the  vertical  offset,  C*,  has  only  a 

slight  effect  on  the  hub  reactions.  Although  not  illustrated  in  figure  20  its  effect  on 
twisting  moment  is  negligible. 

Aerodynamic  load  excitation.  -  The  dynamic  excitation  of  the  rotor  blade  is 
provided  by  a  quasi-steady  representation  of  the  spanwise  airload  distribution  associated 
with  the  cyclic  pitch  environment  of  forward  flight.  As  previously  discussed  the 
aerodynamic  loads  induced  by  the  elastic  deformations  are  also  included.  The 
aerodynamic  load  distributions  consist  of  lift,  drag,  and  pitching  moment.  Values  of  the 
aerodynamic  parameters  used  to  calculate  the  airloads  are  presented  in  Table  3. 

As  previously  discussed  the  unsteady  aerodynamic  forces  which  occur  in  forward 
flight  are  due  to  time  dependent  variations  in  velocity,  angle  of  attack,  and  elastic 
deflections  encountered  by  the  rotating  blade.  The  quasi-steady  representation  provides 
1/Rev,  2/Rev,  and  3/Rev  excitation  airloads.  The  amplitude  of  the  airloads  decreases 
as  the  harmonic  order  increases.  For  each  of  the  three  harmonics,  the  amplitude  of  the 
vertical  force  is  much  larger  than  the  amplitude  of  the  in-plane  force. 

For  airloads  which  are  identical  on  each  rotor  blade,  the  only  loads  and  moments 
which  the  blades  may  transmit  to  the  rotor  hub  are  those  which  are  integral  multiples 
(N)  of  the  number  of  blades  (ref.  17).  The  N/Rev  rotor  hub  loads  are  the  greatest 
contributors  to  helicopter  vibration.  The  source  of  N/Rev  hub  loads  are  N-l,  N,  N+l 
rotor  blade  airload  harmonics.  Therefore  one  solution  for  reducing  N/Rev  vibration  is 
to  cancel  or  eliminate  these  harmonics  of  blade  loading.  As  discussed  above,  the  quasi- 
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TABLE  3.  -  VALUES  OF  AERODYNAMIC  PARAMETERS  FOR  THE 


NONUNIFORM  BLADE  IN  FORWARD  FLIGHT 

Advance  ratio,  y  j .  0.3 

Shaft  angle,  a$,  deg .  -5 

Cosine  cyclic  pitch  coefficient,  6  jc,  deg . -1.5 

Sine  cyclic  pitch  coefficient,  6ls’  deg .  3.0 

Collective  pitch,  0c,  deg .  13 

Pretwist,  0  ,  deg .  * 9 

Airfoil  lift-curve  slope,  a,  per  rad .  2  tt 

Profile  drag  coefficient,  c^ . 0.01 

Precone  angle,  B  ,  rad .  0.1 

pc 

Blade  rotational  speed,  fl  ,  rpm .  289 
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steady  representation  of  the  airloads  contains  only  three  harmonics.  In  the  present 
analysis,  the  effectiveness  of  the  flapping  pendulum  in  reducing  the  blade  root  forces 
due  to  2/Rev  and  3/Rev  airloads  is  discussed.  Throughout  the  following,  unless  stated 
otherwise,  all  numerical  results  pertain  to  the  nonuniform  blade. 

The  root  reactions  of  the  nonuniform  blade  with  a  flapping  pendulum  are  computed 
for  a  2/Rev  spanwise  airload  distribution.  The  results  for  the  vertical  shear  at  the  hub 
are  plotted  in  figure  21.  The  pendulum  weight  is  133.452  N.  From  the  hub  reactions  it 
is  observed  that  the  pendulum  of  uncoupled  natural  frequency  2.02/Rev  located  at  22% 
of  the  span  is  a  proper  design  because  it  attenuates  the  hub  vertical  shear  to  22%  and 
the  out-of-plane  moment  to  31%.  This  choice  of  pendulum  attenuates  the  in-plane 
moment  by  only  4%  and  amplifies  the  in-plane  force  by  66%.  However  for  the  2/Rev 
airload  the  in-plane  hub  force  is  small  compared  to  the  hub  vertical  shear. 

For  the  optimum  pendulum  tuning  (  ui  =  2.02/Rev  and  =  22%  span)  the  root 
reactions  are  computed  for  different  values^ of  the  pendulum  mass  and  the  results  are 
presented  in  figure  22.  It  car  be  observed  that  a  94.53  N  pendulum  weight  (0.71  weight 
ratio)  is  a  proper  choice  beca  e  it  attenuates  the  hub  shear  to  2.9  %.  It  also  attenuates 
the  hub  out-of-plane  moment  to  17%. 

To  establish  the  importance  of  representing  the  excitation  as  a  distributed  loading 
the  responses  due  to  the  distributed  airload  are  compared  with  corresponding  responses 
associated  with  a  concentrated  harmonic  force  applied  at  the  tip  with  the  same 
frequency.  Figure  23  shows  the  results  for  a  2/Rev  concentrated  out-of-plane  force  at 
the  tip.  It  can  be  noticed  from  figure  23  that  the  flapping  pendulum  effect  on  the  hub 
vertical  shear  is  similar  to  the  results  obtained  for  the  2/Rev  spanwise  airload 
distribution  (fig.  21).  For  specific  values  of  the  pendulum  frequency  and  spanwise 
location,  the  pendulum  attenuates  or  amplifies  the  root  vertical  shear  the  same  way  for 
both  the  distributed  airload  and  the  concentrated  load.  It  is  observed  that  an  optimum 
tuning  occurs  at  22%  spanwise  location  with  pendulum  frequency  equal  to  2.02/Rev. 
This  optimum  tuning  attenuates  the  vertical  shear  and  the  out-of-plane  moment  to  17 
and  38%,  respectively.  It  amplifies  the  in-plane  force  by  63%.  It  may  be  recalled  that 
the  optimum  pendulum  for  the  2/Rev  distributed  airloads  also  occurs  at  x.  =  22%  span, 
u)  =  2.02/Rev,  and  W  =  94.53  N.  A 

The  pendulum  mass  optimization  for  the  2/Rev  concentrated  load  is  shown  in  figure 
24.  It  can  be  observed  that  a  100.09  N  pendulum  weight  (0.75  weight  ratio)  is  a  proper 
choice  because  it  attenuates  the  hub  shear  to  8%.  It  also  attenuates  the  root  out-of¬ 
plane  moment  to  26%. 

For  the  2/Rev  concentrated  load  at  the  tip,  discussed  above,  the  elastic- 
aerodynamic  coupling  is  not  included.  However,  it  is  important  to  study  the  effect  of 
these  terms  when  the  blade  is  excited  by  a  concentrated  load.  Figure  25  shows  the 
results  for  a  2/Rev  concentrated  load  at  the  tip  with  the  aerodynamic  coupling  terms 
included,  as  well  as  the  results  without  these  terms  (W  =  133.45  N).  It  can  be  noted 
from  figure  25  that  for  both  cases  the  root  vertical  shear  has  the  same  characteristics. 
However,  the  amount  of  attenuation  is  different.  Figure  25  shows  that,  for  the  case 
with  aerodynamic  coupling  the  optimum  tuning  occurs  at  x^  =  26%  span  and 
<*>  -  2.02/Rev.  This  choice  of  pendulum  attenuates  the  root  out-of-plane  shear  and 

moment  to  38%  and  84%,  respectively.  The  results  of  varying  the  pendulum  mass 
indicate  that  the  reaction  attenuation  varies  only  by  1.0%.  The  results  also  show  that  a 
1  33.45  N  is  a  proper  choice.  Comparing  the  pendulum  tuning  for  the  2/Rev  distributed 
airloads,  the  2/Rev  concentrated  load,  and  the  2/Rev  concentrated  load  with 
aerodynamic  coupling,  it  can  be  concluded  that  a  pendulum  tuned  for  a  concentrated 
simple  harmonic  load  acting  at  the  blade  tip  can  serve  the  purpose  of  a  distributed 
simple  harmonic  loading  at  the  same  frequency  as  far  as  the  pendulum  tuning  and  mass 
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optimization  is  concerned.  However,  it  is  necessary  to  utilize  appropriate  distributed 
airloads  to  accurately  determine  the  attenuation  of  the  root  reactions. 

Figure  26  shows  the  effect  of  pendulum  tuning  on  the  root  vertical  shear  for  3/Rev 
spanwise  airload  distribution.  The  results  of  figure  26  indicate  that  a  flapping  pendulum 
of  frequency  3.10/Rev  located  at  26%  of  the  span  is  an  optimum  tuning.  This  choice  of 
pendulum  attenuates  the  root  vertical  shear  and  the  out-of-plane  moment  to  4.2  and 
26.8%,  respectively.  It  attenuates  the  root  in-plane  moment  only  to  97%.  However,  the 
in-plane  force  at  the  hub  is  amplified  by  44%.  This  adverse  effect  is  the  penalty  being 
paid  for  achieving  the  drastic  reductions  in  the  root  out-of-plane  force  and  moment. 

The  root  reactions  are  computed  for  various  masses  of  the  optimum  pendulum. 
Results  for  the  vertical  shear  are  presented  in  figure  27.  These  data  show  that  a  155.69 
N  pendulum  weight  (1.17  weight  ratio)  is  a  proper  choice.  For  this  weight  the  vertical 
shear  and  out-of-plane  moment  are  attenuated  to  3.3  and  25.2%,  respectively. 

The  above  results  for  a  simple  flapping  pendulum  on  the  nonuniform  blade  are 
summarized  in  Table  4. 


Simple  Lead-Lag  Pendulum 

As  discussed  in  the  previous  section,  the  simple  flapping  pendulum  can  be  tuned  to 
achieve  significant  reductions  in  the  root  out-of-plane  forces.  To  significantly 
attenuate  the  in-plane  forces  a  simple  lead-lag  pendulum  is  needed.  In  this  section,  the 
effect  of  mounting  a  lead-lag  pendulum  on  a  rotor  blade  is  discussed.  The  excitation 
loads  to  be  considered  are  4/Rev  concentrated  load  at  the  blade  tip  and  2/Rev  and 
3/Rev  spanwise  distributed  airloads.  Throughout  the  following  all  numerical  results 
pertain  to  the  nonuniform  blade. 

Concentrated  load  excitation.  -  The  root  reactions  of  the  nonuniform  blade  with  a 
simple  lead-lag  pendulum  are  computed.  The  results  for  the  in-plane  force  amplitude  at 
the  hub  are  plotted  in  figure  28.  The  study  of  the  effect  of  the  pendulum  tuning  on  the 
root  in-plane  reaction  is  of  prime  importance  since  the  lead-lag  pendulum  attenuates 
mainly  the  in-plane  forces.  The  pendulum  weight  is  equal  to  133.4  52  N.  From  the 
results  it  can  be  observed  that  the  pendulum  of  uncoupled  natural  frequency  3.85/Rev 
located  at  26%  ot  the  span  is  a  proper  design  because  it  attenuates  the  hub  in-plane 
force  to  5%  and  the  in-plane  moment  to  44.6%.  This  choice  of  pendulum  also 
attenuates  the  root  vertical  shear  and  the  out-of-plane  moment  to  20.6  and  16.8%, 
respectively.  In  principle,  it  is  feasible  to  design  a  lead-lag  simple  pendulum  to  achieve 
significant  reductions  in  the  hub  in-plane  force  and  moment  without  producing  an 
adverse  effect  on  the  out-of-plane  forces. 

For  the  optimum  lead-lag  pendulum  (ui  =  3.85/Rev  and  =  26%  span)  the  root 
forces  are  computed  for  different  values  of  tfte  pendulum  mass.  The  results  are  shown 
in  figure  29.  It  can  be  observed  from  figure  29  that  if  the  pendulum  weight  ratio  is 
reduced  from  1.0  to  0.3,  then  the  normalized  in-plane  force  at  the  hub  increases  from 
5%  to  39%.  This  result  indicates  that  even  though  the  uncoupled  natural  frequency  of 
the  lead-lag  pendulum  is  independent  of  the  mass  of  the  pendulum,  a  proper  choice  of 
the  mass  is  required  to  minimize  the  root  forces.  It  may  be  recalled  that  this 
conclusion  was  also  observed  for  the  simple  flapping  pendulum. 

It  is  appropriate,  at  this  stage,  to  compare  the  results  for  the  optimum  tuned 
flapping  and  lead-lag  pendulums.  Shown  in  Table  5  are  the  optimum  pendulum 
parameters  for  both  the  flap  and  lead-lag  pendulum  configurations.  Also  shown  in  this 
table  are  the  normalized  root  forces  of  the  nonuniforrn  blade  due  to  a  4/Rev 
concentrated  load  at  the  tip.  From  Table  5  it  can  be  observed  that  for  both  the  flap  and 
lead-lag  pendulums  the  uncoupled  frequencies  are  the  same.  Also  the  uncoupled 
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TABLE  5.  -  ROOT  REACTIONS  FOR  OPTIMUM  FLAP  AND 
LEAD  -  LAG  PENDULUMS  AT  4/REV 


Normalized  root  reaction, 

Flapping 

Lead  -  lag 

percent 

pendulum3 

pendulum** 

M 

139.0 

85.0 

X 

M 

5.8 

16.8 

y 

M 

172.0 

44.6 

z 

V 

43.6 

5.0 

y 

V 

0.2 

20.6 

z 

a  a.  =  3.85/Rev  }  xA 
P  A 

=  22%  span  ; 

W  =  66.72  N 

**  u»  =  3.85/Rev  ;  x, 

P  A 

=  26%  span  ; 

W  =  133.45  N 

L 
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frequency  (3.85/Rev)  is  in  the  vicinity  of  the  external  load  frequency  (4/Rev).  It  may  be 
concluded  that,  in  general,  the  optimum  pendulum  frequency  depends  mainly  or.  the 
excitation  frequency.  This  observation  is  due  to  the  fact  that  the  optimum  tuned 
pendulum  attenuates  the  hub  reactions  by  eliminating  the  resonant  responses  of  the 
blade  (in  addition  to  generating  an  appropriate  reaction  force  at  the  hinge).  The 
optimum  pendu<"m  eliminates  these  resonant  responses  by  displacing  the  natural 
frequencies  tha  e  in  the  vicinity  of  the  excitation  frequency,  and  since  the  pendulum 
is  more  effect  .  in  altering  the  blade  frequencies  adjacent  to  the  uncoupled  pendulum 
frequency  (ref.  7),  then  the  uncoupled  pendulum  frequency  must  be  close  to  the 
excitation  frequency.  As  expected,  the  flapping  pendulum  significantly  attenuates  the 
hub  vertical  shear,  whereas  the  lead-lag  pendulum  reduces  the  hub  in-plane  force 
drastically  as  shown  in  Table  5.  Finally,  it  can  be  observed  from  Table  5  that  for  this 
particular  blade  the  flapping  pendulum  produces  an  adverse  effect  on  the  twisting  and 
in-plane  moments,  whereas  the  lead-lag  pendulum  attenuates  the  twisting  and  out-of¬ 
plane  moments. 

The  spanwise  variations  of  deflections  and  forces  for  the  nonuniform  blade 
corresponding  to  the  optimum  lead-lag  pendulum,  which  is  located  at  26%  span,  are 
shown  in  figure  30  together  with  the  values  obtained  without  the  pendulum.  The  lead- 
lag  and  flapwise  deflections  are  normalized  with  respect  to  the  chord  (c  =  0.53  m).  Only 
the  magnitude  of  the  forces  and  moments  are  plotted.  It  can  be  observed  from  figure 
30  that  the  lead-lag  pendulum  redistributes  the  structural  loads  of  the  blade  such  that 
the  hub  forces  are  attenuated.  The  pendulum  decreases  the  loads  inboard  of  its  location 
and  increases  the  loads  outboard.  Figure  30  shows  that  the  lead-lag  deflection  of  the 
blade-pendulum  system  is  significantly  reduced.  Also,  it  can  be  observed  that  the  in¬ 
plane  force,  V  ,  is  characterized  by  a  large  change  at  the  lead-lag  pendulum  attachment 
point  due  to  t¥ie  internal  reaction  of  the  lead-lag  pendulum  mass.  It  may  be  recalled 
that  the  flapping  pendulum  produces  a  large  change  in  the  vertical  shear.  These 
reactions  at  the  attachment  point  redistribute  the  loads  on  the  blades  in  such  a  way  that 
only  a  small  portion  of  the  shears  and  moments  are  transmitted  to  the  hub. 

A  parametric  study  is  conducted  for  the  optimum  lead-lag  pendulum 
(ui  =  3.85/Rev  and  x  .  =  26%  span).  The  parameters  varied  include  prepitch,  pretwist, 
pre^cone,  and  pendulum  ninge  offsets. 

Figure  31  shows  the  effect  of  the  root  collective  pijch  (prepitch)  on  the  hub 
reactions.  The  pretwist  angle  for  these  results  is  equal  to  -9  .  It  can  be  observed  from 
figure  31  that  increasing  the  prepitch  decreases  the  hub  m-plane  force  to  a  minimum 
value  after  which  it  increases.  The  in-plane  hub  moment  is  increased  by  increasing  the 
prepitch  (in  the  normal  range).  Figure  31  shows  that  the  prepitch  has  a  negligible  effect 
on  both  the  out-of-plane  shear  and  moment  at  the  hub. 

The  effect  of  pretwist  on  the  hub  reactions  is  shown  in  figure  32.  The  root 
collective  pitch  for  these  results  is  equal  to  13°.  Figure  32  shows  that  both  the  hub  in¬ 
plane  force  and  moment  are  increased  by  decreasing  the  blade  pretwist.  Thus,  the  lead- 
lag  pendulum  effectiveness  in  attenuating  the  root  in-plane  force  and  moment  is 
reduced  by  decreasing  the  blade  pretwist.  It  can  be  observed  from  figure  32  that  the 
effect  of  the  pretwist  on  both  the  out-of-plane  force  and  moment  at  the  hub  is  small. 

The  root  reactions  as  a  function  of  the  blade  precone  angle  are  computed,  and  the 
results  are  presented  in  figure  33.  It  can  be  observed  that  the  precone  angle  has  no 
significant  influence  on  the  hub  reactions.  It  may  be  recalled  that  the  same  conclusion 
was  observed  for  the  flapping  pendulum. 

The  results  from  varying  the  lead-lag  pendulum  hinge  chordwise  offset  (n  *)  are 
shown  in  figure  34.  It  may  be  noted  that  the  optimum  lead-lag  pendulum  is  tuned  with 
its  hinge  located  at  the  shear  center  of  the  blade  cross  section  (  n  A  =  0,  =  0).  It  can 
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be  observed  from  figure  34  that  the  pendulum  hinge  chordwise  offset  has  only  a  slight 
effect  on  both  the  hub  in-plane  force  and  moment.  However,  both  the  hub  out-of-plane 
force  and  moment  are  decreased  sigificantly  to  practically  zero,  after  which  they  are 
increased,  by  increasing  the  chordwise  offset. 

Figure  35  shows  the  root  reactions  for  different  values  of  the  lead-lag  pendulum 
hinge  vertical  offset,  £  ..  This  offset  has  a  small  effect  on  the  root  in-plane  moment. 
However,  the  in-plane  force  at  the  hub  is  increased  by  increasing  the  vertical  offset  as 
shown  in  figure  35.  Also,  both  the  root  out-of-plane  force  and  moment  are  increased  by 
increasing  the  lead-lag  pendulum  hinge  offset,  £  It  can  also  be  noted  that  both  the 
chordwise  and  vertical  offsets  have  a  negligible  effect  on  the  hub  twisting  moment 
reaction. 

Aerodynamic  load  excitation.  -  The  root  reactions  of  the  nonuniform  blade  with  a 
lead-lag  pendulum  are  computed  for  a  3/Rev  spanwise  airload  distribution.  The  results 
for  the  root  in-plane  force  are  shown  in  figure  36.  It  can  be  observed  from  figure  36 
that  for  this  particular  load  the  lead-lag  pendulum  attenuates  the  root  in-plane  force 
only  by  a  small  amount.  However,  the  root  in-plane  force  is  small  for  all  three 
harmonics  of  the  spanwise  airload  distribution.  The  lead-lag  pendulum  of  uncoupled 
frequency  3.10/Rev  located  at  19%  of  the  span  is  a  proper  choice  because  it  attenuates 
the  hub  in-plane  force  to  65%.  However,  it  amplifies  the  other  root  reactions  by  5  to 
25%. 

The  effect  of  lead-lag  pendulum  tuning  on  the  root  in-plane  force  for  a  2/Rev 
spanwise  airload  distribution  and  a  2/Rev  concentrated  out-of-plane  load  is  shown  in 
figures  37  and  38.  It  can  be  observed  from  these  figures  that  for  specific  values  of  the 
pendulum  frequency  and  spanwise  location,  the  pendulum  attenuates  or  amplifies  the 
root  in-plane  force  the  same  way  for  both  the  distributed  airload  and  the  concentrated 
load.  The  pendulum  of  frequency  2.01/Rev  located  at  9%  span  attenuates  the  root  in¬ 
plane  force  to  60%  for  the  distributed  airload  and  42%  for  the  concentrated  load  at  the 
tip.  As  concluded  earlier,  for  the  flapping  pendulum,  a  pendulum  can  be  tuned  for  a 
concentrated  simple  harmonic  load  at  the  tip  as  well  as  by  a  distributed  simple 
harmonic  loading  of  the  same  frequency. 

CONCLUDING  REMARKS 

Comprehensive  analytical  design  procedures  for  the  installation  of  simple 
pendulums  on  the  blades  of  a  helicopter  rotor  to  suppress  the  root  reactions  were 
developed.  To  achieve  this  goal,  a  frequency  response  analysis  of  typical  rotor  blades 
excited  by  a  harmonic  variation  of  spanwise  airload  distributions  as  well  as  a 
concentrated  load  at  the  tip  was  conducted.  A  single  nonuniform  rotor  blade  with 
hingeless  hub  restraint  undergoing  coupled  flapwise  bending,  chordwise  bending,  and 
torsional  vibrations  was  considered.  The  equa*  ions  of  motion  of  both  flap  and  lead-lag 
pendulums  were  derived.  Expressions  for  the  aerodynamic  loads  associated  with  the 
elastic  deformations  and  the  cyclic  pitch  environment  of  forward  flight  were  presented. 
The  transfer  matrix  method  was  used  to  determine  the  optimum  pendulum  tuning  and 
mass  to  suppress  the  hub  reactions.  Lastly,  the  effects  of  various  structural  dynamic 
properties  of  the  blade  on  the  optimum  pendulum  configuration  were  established. 

The  results  of  the  study  indicate  that  the  following  conclusions  can  be  di  awn. 

(1)  In  general,  it  is  feasible  to  design  a  simple  pendulum  to  suppress  the  vibratory 
loads  at  the  hub  of  a  rotor  blade.  Further: 

(i)  A  properly  designed  flapping  plane  simple  pendulum  attenuates  the  root  out- 
of-plane  force  and  moment  significantly.  In  principle,  the  reductions  in  the  root  out-of¬ 
plane  reactions  can  be  achieved  without  producing  adverse  effects  on  the  root  in-plane 
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forces.  However,  in  some  of  the  cases  investigated,  the  hub  in-plane  forces  are 
amplified  with  the  flapping  pendulum.  This  is  a  design  penalty  which  results  from 
reducing  the  out-of-plane  reactions. 

(ii)  An  optimum  lead-lag  pendulum  design  attenuates  the  root  in-plane  force  and 
moment. 

(iii)  For  the  cases  investigated,  the  flapping  plane  simple  pendulum  is  more 
effective  in  attenuating  the  root  out-of-plane  reactions  than  the  lead-lag  pendulum  in 
suppressing  the  in-plane  forces  at  the  hub. 

(2)  A  properly  tuned  pendulum  can  attenuate  the  vibratory  loads  in  two  ways:  (1) 
by  eliminating  the  resonant  responses  of  the  blade  by  displacing  the  initial  frequencies 
that  are  in  the  vicinity  of  the  excitation  frequency,  and  (2)  by  generating  appropriate 
forces  at  its  attachment  point  with  the  blade.  These  forces  redistribute  the  loads  on 
the  blade  so  that  only  a  small  portion  of  the  shears  and  moments  are  transmitted  to  the 
hub. 

(3)  For  optimum  tuning,  the  parameters  to  be  determined  are  the  pendulum 
uncoupled  natural  frequency,  the  pendulum  spanwise  location,  and  its  mass.  Further: 

(i)  The  pendulum  uncoupled  natural  frequency  is  a  function  of  the  (a)  rotational 
speed  of  the  blade,  (b)  spanwise  location  of  the  pendulum,  (c)  its  length,  (d)  steady-state 
deflection  of  the  pendulum,  and  (e)  blade  pitch  angle  at  the  pendulum  hinge  location. 

(ii)  The  optimum  pendulum  frequency  is  in  the  vicinity  of  the  excitation 
frequency.  In  general,  both  the  flapping  plane  pendulum  and  the  lead-lag  pendulum  have 
the  same  optimum  frequency  for  the  same  excitation  frequency.  In  the  case  of  an  off- 
tuned  pendulum  the  results  range  from  a  slight-to-moderate  degradation  in  pendulum 
absorber  effectiveness  to  a  drastic  increase  in  hub  vibration. 

(iii)  The  uncoupled  natural  frequency  of  the  pendulum  is  independent  of  the 
mass  of  the  pendulum,  and  a  proper  choice  of  the  mass  is  required  to  generate  an 
optimum  force  to  minimize  the  root  reactions. 

(4)  A  pendulum  can  be  tuned  and  its  optimum  mass  determined  by  excitation  with  a 
concentrated  simple  harmonic  load  at  the  tip.  However,  it  is  necessary  to  utilize 
appropriate  distributed  airloads  to  accurately  determine  the  attenuation  of  the  root 
reactions. 

(5)  Damping  at  the  pendulum  hinge  has  a  very  small  effect  on  the  hub  reactions 
once  the  optimum  pendulum  tuning  is  established. 

(6)  The  parametric  study  for  a  4/Rev  flapping  plane  pendulum  and  a  4/Rev  lead-lag 
pendulum  indicates: 

(i)  Blade  precone  has  no  significant  effect  on  the  hub  reactions  for  both  types  of 
pendulum. 

(ii)  The  effect  of  collective  pitch  on  the  hub  vertical  shear  is  small,  but  has  an 
appreciable  effect  on  the  hub  out-of-plane  moment  for  a  flapping  plane  pendulum.  For 
the  lead-lag  pendulum  the  collective  pitch  has  a  negligible  effect  on  both  the  hub  out- 
of-plane  force  and  moment.  Increasing  the  collective  pitch,  the  hub  in-plane  force  and 
moment  are  significantly  raised  in  the  case  of  the  flapping  pendulum,  whereas  they  are 
lowered  to  a  minimum  value  after  which  they  are  increased  in  the  case  of  a  lead-lag 
pendulum. 

(iii)  The  pendulum  effectiveness  in  attenuating  the  hub  in-plane  force  and 
moment  is  reduced  by  decreasing  the  blade  pretwist  for  both  types  of  pendulum.  In  the 
case  of  the  flapping  pendulum,  pretwist  has  a  significant  effect  on  the  hub  out-of-plane 
moment. 

(iv)  The  pendulum  hinge  chordwise  offset  has  a  small  effect  on  the  hub  out-of¬ 
plane  reactions  in  the  case  of  a  flapping  plane  pendulum.  Also,  for  this  type  of 
pendulum  the  hub  in-plane  reactions  are  slightly  decreased  by  increasing  the  chordwise 
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offset.  In  the  case  of  a  lead-lag  pendulum,  the  hinge  chordwise  offset  has  only  a  slight 
effect  on  the  root  in-plane  reactions,  whereas  the  out-of-plane  forces  are  significantly 
decreased  by  increasing  the  chordwise  offset. 

(v)  The  pendulum  hinge  offset  normal  to  the  chord  has  only  a  slight  effect  on 
the  hub  reactions  for  the  flapping  plane  pendulum.  In  the  case  of  a  lead-lag  pendulum, 
the  hub  in-plane  force  and  both  the  out-of-plane  force  and  moment  are  increased  by 
increasing  the  hinge  offset.  But  the  hinge  offset  has  only  a  small  effect  on  the  hub  in¬ 
plane  moment. 

(vi)  For  both  the  flapping  and  lead-lag  pendulum  configurations  the  pendulum 
hinge  offset  in  both  the  chordwise  and  normal  directions  had  a  negligible  effect  on  the 
twisting  moment  reaction  at  the  hub. 
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APPENDIX  A 


COORDINATE  SYSTEMS  AND  TRANSFORMATIONS 


Several  coordinate  systems  are  used  in  deriving  the  blade  and  pendulum  equations 
of  motion.  Those  which  are  common  to  both  the  dynamic  and  aerodynamic  aspects  of 
the  derivation  are  shown  in  figures  A-l  and  A-2.  The  orthogonal  axes  system  X,  Y,  Z 
(fig.  A-l)  and  associated  unit  vectors  T,  J,  K  are  fixed  in  an  inertial  frame  .  It  is 
assumed  that  the  rotor  hub  has  neither  translational  nor  angular  motion.  Orthogonal 
axes  x  ,  y,  Z  are  fixed  in  a  ^ference  frame  QS  which  rotates  with  respect  to  Qft  at 
constaRt  angular  velocity  ftK.  Point  O,  which  is  common  to  both  the  inertial  and 
rotating  frames,  is  located  at  the  root  of  the  beam.  The  x  axis,  which  lies  along  the 
elastic  axis  of  the  undeformed  beam,  is  inclined  to  the  plane  of  rotation  (and  to  the  x 
axis)  by  ^he  jjrecone  angle  8  .  The  orthogonal  axes  x,  y,  z  and  the  corresponding  uniR 

vectors  1,  j,  k,  therefore,  arePalso  fixed  in  .  Beam  bending  deformations  shown  in 
figure  A-2  for  9=0  are  described  by  the  displacements  of  the  elastic  axis  u,  v,  w 
parallel  to  1,  j,  &  respectively. 

The  transformation  between  the  X,  Y,  Z  coordinates  of  the  inertial  reference 
frame,^,  and  the  x  ,  y,  Z  coordinates  of  the  rotating  reference  frame, is  given  by 


X 

Y 


cos  ft  t 
sin  ft  t 
0 


(A-l  ) 


Rotation  of  the  x  ,  y,  Z  coordinate  system  by  the  precone  angle  P  about  the  y  axis 
yields  the  x,  y,  z  Rindeformed  coordinate  system.  These  two  systemsPare  related  by  the 
following  transformation. 


°os  B^  0 


0 

sin  B 


pc 


1 

0 


-sin  B 
0 

cos  B 


pc 


pc 


(A-2) 


For  small  values  of  the  precone  angle  equation  (A-2)  can  be  approximated  as 


(A-  3 ) 


The  relationship  between  the  deformed  x’,  y’,  z'  coordinates  and  the  undeformed  x, 
y,  z  coordinates  illustrated  in  figure  A-2  has  been  developed  in  refer^n^e  10.  A 
transfor  ium  matrix,  [  T  ]  ,  which  relates  their  respective  unit  vectors  i',  j', Ur' and  i, 
I,  £,  can  .  defined  as 
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APPENDIX  B 


ACCELERATION  COMPONENTS 

The  purpose  of  this  appendix  is  to  formulate  the  components  of  absolute 
acceleration  as  measured  in  the  deformed  coordinate  system  for  a  general  point 
(x.,y.,z.)  in  the  undeformed  system.  It  is  desired  to  use  the  deformed  coordinate 
system  since  it  is  easier  to  formulate  the  pendulum  equilibrium  equations  in  this  system. 
These  acceleration  components  are  to  be  described  in  terms  of  section  coordinates  and 
the  elastic  blade  displacements  and  rotations  to  the  first  order. 

A  position  vector  in  the  inertial  reference  system  can  be  written  as 


xl 


Y5 


2  £ 


(B-l  ) 


and  the  corresponding  acceleration  components  are 


X  ;  a, 


Y  ;  a- 


(B-2) 


The  components  of  acceleration  in  the  rotating,  nonpreconed  system  can  be  obtained 
from  the  component  transformations  in  equation  (A-l)  as 


xp 


X  cos  ft  t  +  Y  sin  ft  t 

Y  cos  ft  t  -  X  sin  ft  t 

"z 


(B-3) 


Substitution  of  equation  (B-3)  into  equation  (A-3)  yields  the  acceleration  components  in 
the  undeformed  system  which  is  rotating  at  a  small  precone  angle  as  illustrated  in 
figure  A-l.  These  components  become 


X  cos  ftt  ♦  Y  sinftt  +  Z  3 

Y  cos  ftt  -  X  sin  ft  t 

Z  -  (X  cos  ft  t  +  Y  sin  ftt)  g 

PC 


pc 


(B-4 ) 


The  position  vector  of  equation  (B-l)  can  be  written  for  a  general  point  (x , ,  y.f 
Zj)  in  the  undeformed  system  by  successive  applications  of  equations  (A-l)  and  (A-3)  as 

r  =  (x^  cos  fit  -yj  sin  ftt  -Zj  g  cos  ftt)T 

♦(x.  sin  ftt  +y.cos  ftt  -z,  g  sin  ft  t)  3 
1  '1  1  pc 


h(xi  v  +  2i>* 


(B-5) 
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The  second  time  derivative  of  the  components  of  equation  (B-5)  are 

X  =  x  j  cos  ft  t  ~y  j  sin  ft  t  -z  j  0  pc  cos  ft  t  ( B  -  6  ) 

-2  ft  x  i  sin  ftt  -2  ft  y  j  cos  ft  t  +  2  ft  z  ^  0  pc  sin  ft  t 

-  ft2  x .  cos  ft  t  +  ft2  y .  sin  ft  t  +  ft2  z .  0  cos  ft  t 

I  1  1  pc 

Y  =  Xj  sin  ftt  +■  y^  cos  ft  t  -z^  0pc  sin  ftt  (B-7) 

+  2  ft  x .  cos  ft  t  -2  ft  y .  sin  ft  t  -2  ft  z .  0  cos  ft  t 

I  1  1  pc 

2  2  2 
-ft  Xj  sin  ft  t  -  ft  y^  cos  ft  t  +  ft  zi  0pc  sin  ft  t 

Z  =  Xi  0pc  +  z.j  (B-  8) 


Substituting  these  derivatives  into  the  acceleration  components  of  equation  (B-4)  and 
neglecting  second  degree  terms  in  the  precone  angle  yields 

ax  =  xt  -2ft  yt  -  ft2xj  ♦  ft2  z,  0pc  j 

ay  =  y!+2ftXj  -2ftZj  0pc  -  ft2  y}  \  ( B-9 ) 

az  =  Z1  +2nyl  Bpc  +  X1  V  J 

In  the  deformed  body,  the  coordinates  Xj,yj,Zj  of  the  general  point  in  the 
undeformed  coordinate  system  can  be  expressed  in  terms1  of  the  elastic  deflections  u,  v, 
w,  <J>  and  the  section  coordinates  n,C  (same  as  x',  y',  z')  with  the  aid  of  equation  (A- 
4)  as 


-  x{  - 

yl 

= 

■  X  ♦  u" 

V 

*  [T]T 

1 

- 1 

-Z1  - 

w 

These  coordinates  become 


Xj  -  x 


+  u  ♦  £  -v'  cos  (0  +  $  )  -  £  sin  (  6  +  ♦  >] 

-w'  £q  sin  (0  +  <f>)  +  £  cos  (0  +  <J>  )  J 


( B-  1 0 ) 


(B-ll  ) 
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yt  =  V  +  v'  £  +  n  COS  (0  +  4) )  -  c  sin  (e  +  4>  )  (B-12) 

Z  J  =  w  ♦  w*  £  +  n  Sin  ( e  +  $  )  +  c  cos  (6  +  <f> )  ( B-  1 3  ) 


When  these  coordinates  and  their  time  derivatives  are  substituted  into  equation  (B-9) 
the  acceleration  components  to  first  order  become 


a 


x 


a 


y 


a 


z 


£  -2  fi  v  -  2  (  x  +  £  ) 

+  2  ;  sin  (0  +  <$>  )  -2  ft  q  cos  (  0  +  0  )  (B-  14  ) 

v+2ni-n2v 

2  2 

+  (  ft  S  -  c  )  sin  ( 6  +  «J) )  -  (  Q  n  -  n  )  cos  (e  +  4,  ) 

( B- 1 5 ) 

w+  ^  ^  ®pc  (  x  +  Z  ) 

♦  n  sin  (  0  +  4> )  +  C  cos  (  6  +  )  ( B-  1 6 ) 


In  the  above  acceleration  expressions,  the  coordinates  n  ,C  are  of  a  general  point  off 
the  blade.  These  acceleration  components  in  the  undeformed  system  can  be 
transformed  into  the  defon.-'sd  system  components  by  again  applying  the  transformation 
of  equation  (A-4)  as 


a  ,  "I 
x’ 

a  , 

y 

II 

t - 1 

H 

LJ 

'  a  " 

X 

a 

y 

-  az'  - 

a 

zJ 

(B-17 ) 


The  final  acceleration  components  in  the  deformed  system  can  be  further  simplified  by 
applying  the  trigonometric  approximations  of  equation  (8)  to  yield 


ax,  =  t  -2  ft  v  -  ft  z  (  x  ♦  £)  ♦  2  fl  c  sin  6 


-2  n  n  cos  e 


( b- 1 8 ) 
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APPENDIX  C 


AERODYNAMIC  COUPLING  COEFFICIENTS 

In  equation  (132)  the  elastic  deformations-dependent  aerodynamic  loads  are 
expressed  in  a  matrix  form.  The  nonzero  elements  of  these  matrices  are 


B65 

= 

1  ,Cv 2  0  2 

2  pac  (jj)  ft  x 

B 1 0,3 

= 

1  2 
pac  !1  xR\ 

B10,5 
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-  A  pac  (i  ft  2  x  p.  R  0.  +  ft2  p2  R2  6  ) 

2  2  rf  lc  Hf  pc 

Bli,l 

= 

-  4  p  ac  (-  ft2  X  B  ) 

B1 1,3 

n 

1  ,  n2  2  1  2  2  d2v 

-j  pac(-il  x  -^Pfft  R) 

B1 1,4 

r 

pac(-ft2x  pf  R  0lc-i  ft2  p2R2  Bpc) 

B1 1,5 

= 

-  4  P  ac  (-  |  ft2  x) 

C63 

- 

^  pac  (|)2  (2  ftx) 

C 1 0,2 

-  4  pac  ftx  (  6c  +  0t)  +  |  pf  ftR  6  ls  -2ftR  \ 

Cll,l 

= 

pac  j^2  ft  x  (  0  c  +  0t)-pfftR0ls+ftRx] 

c 

^11,2 

- 

-  j  p  ac  (-ft  x) 

C1 1,3 

-j  pac  (^  ft  x) 
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B  (  = 

si  1,5 

Csl0,2 

C 

si  1,1 

For  n  =  3  ; 

Bc10,5 
Bc  1 1,4 
Bs10,5 
Bs11,4 


Pac  (  -  ^  11 2  UfR2) 
pac  (j  yf  ft  R  9lc) 

-  j  pac  (-  p{  HR  eJc) 

pac(-i  n2  p2r2  els) 
-1  Pac(^2Mf2R2  ei5) 
-i  pac(i  n2  y 2  R2  0lc) 
pac(-i  122y2R2  e,c) 


APPENDIX  D 


ROTOR  CONFIGURATIONS 

Two  rotor  configurations  have  been  considered  in  this  investigation.  The  first  is  a 
uniform  blade,  and  the  second  is  typical  of  an  operational  blade  with  nonuniform 
properties. 

The  uniform  blade  has  the  following  properties: 

Radius,  in . 260 

Linear  twist,  deg .  -10 

Weight  per  unit  length,  lb/in . 0 . 5796 

Flapwise  mass  moment  of  inertia  per  unit  length,  lb-in . 0. 3476 

Chordwise  mass  moment  of  inertia  per  unit  length,  lb-in . 15. 456 

2  8 

Flapwise  bending  stiffness,  lb-in  . 0 . 3x10 

Chordwise  bending  stiffness  lb-in^ . 0 . 1x10^ 

2  8 

Torsional  rigidity,  lb-in . 0.2x10 

Mass  axis  offset,  in .  -0.6 

Rotational  speed,  rpm .  360 


The  nonuniform  blade  has  the  properties  listed  in  the  following  table  and 
illustrated  in  figures  D-l  through  D-6. 


Radius,  in .  288 

Root  cutout,  in .  39 

Precone,  rad .  0.1 

Rotational  speed,  rpm .  289 
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Figure  D-6.  -  Rotor  blade  lead-lag  moment  of  inertia. 
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(b)  Moments. 


Figure  2.  -  Equilibrium  forces  and  moments. 
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Figure  7.  -  Rotor  blade  airfoil  inflow  geometry. 
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Pendulum  spanwise  location/span 


Figure  10  -  Root  vertical  shear  for  4/Rev  concentrated  excitation 
uniform  blade,  W  ,  =  66.73  N,  nA  =  CA  =  0. 


(b)  Flapwise  deflection. 


Figure  12.  -  Deflections  and  loads  for  4/Rev  concentrated  excitation, 
uniform  blade,  x./R  =  .18,  u>  =  3.S7/Rev. 


Figure  12.  -  Continued. 
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Figure  13.  -  Root  vertical  shear  for  4/Rev  concentrated  excitation, 
nonuniform  blade,  W  =  133.43  N,  n.  =  C.  =  0. 


Normalized  root  vertical  shear 


Figure  14.  -  Root  vertical  shear  for  4/Rev  concentrated 
excitation,  nonuniform  blade, 

Wref  =  133.45  N,  nA=  CA  =  0. 


Normalized 


Normalized 


Normalized  root  in-plane  moment,  M 


Figure  17.  -  Pretwist  effect  for  4/Rev  concentrated 
excitation,  nonuniform  blade. 
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Figure  18.  -  Precone  effect  for  4/Rev  concentrated 
excitation,  nonuniform  blade. 
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Figure  24.  -  Root  vertical  shear  for  2/Rev  concentrated 
excitation,  nonuniform  blade, 


Normalized  root  vertical  shear 


Pendulum  spanwise  location/span 


Figure  25.  -  Root  vertical  shear  for  2/Rev  concentrated  excitation, 
with  and  without  aerodynamic  coupling,  nonuniform 
blade,  W=  133.45  N,  n  A  =  r  A  =  0. 


Normalized  root  vertical  shear 


Figure  26.  -  Root  vertical  shear  for  3/Rev  distributed  airload  excitation, 
nonuniform  blade,  W  =  1 33.45  N,  =  0. 
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Flapping  Pendulum 


Normalized  root  in-ptane  shear 


Figure  28.  -  Root  in-plane  shear  for  4/Rev  concentrated  excitation, 
nonuniform  blade,  W  =  133.45  N,  ri  ^  =  5  ^  =  0. 
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Figure  30.  -  Deflections  and  loads  for  4/Rev  concentrated  excitation, 
nonuniform  blade,  x./R  =  .26,  u>  =  3.85/Rev. 


(c)  Torsional  deflection. 
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(d)  Twisting  moment  (magnitude). 


Figure  30.  -  Continued. 


Figure  30.  -  Continued. 
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Figure  30.  -  Concluded. 
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Figure  35.  -  Pendulum  vertical  hinge  offset,  £  . ,  effect  for  4/Rev 
concentrated  excitation,  nonuniform  blade. 
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Normalized  root  in-plane  shear, 


Normalized  root  in-plane  shear 
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Lead-Lag  Pendulum 


Pendulum  spanwise  location/span 


Figure  36.  -  Root  in-plane  shear  lor  3/Rev  distributed  airload  excitation, 
nonuniform  blade,  W  =  133.45  N,  =  C^?0. 
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Normalized  root  in-plane  shear 


Figure  37.  -  Root  in-plane  shear  for  2/Rev  distributed  airload  excitation, 

*A  =  c  A 


nonuniform  blade,  W  =  133.45  N,  =  C  »  =  0. 
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